arXiv:1508.02569vl [gr-qc] 11 Aug 2015 


The theory of Hawking radiation in laboratory 

analogues 

Scott J Robertson 
Dipartimento di Fisica, 

Universita degli Studi di Pavia, Via Bassi 6, 27100 Pavia, Italy 
email: scott.robertson@unipv.it 

Monday 2 nd July, 2012 


Abstract 

Hawking radiation, despite being known to theoretical physics for nearly 
forty years, remains elusive and undetected. It also suffers, in its original 
context of gravitational black holes, from practical and conceptual difficulties. 
Of particular note is the trans-Planckian problem, which is concerned with the 
apparent origin of the radiation in absurdly high frequencies. In order to gain 
better theoretical understanding and, it is hoped, experimental verification of 
Hawking radiation, much study is being devoted to laboratory systems which 
use moving media to model the spacetime geometry of black holes, and which, 
by analogy, are also thought to emit Hawking radiation. These analogue 
systems typically exhibit dispersion, which regularizes the wave behaviour at 
the horizon at the cost of a more complicated theoretical framework. This 
tutorial serves as an introduction to Hawking radiation and its analogues, 
developing the moving medium analogy for black holes and demonstrating 
how dispersion can be incorporated into this generalized framework. 
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1 Introduction 


Black holes - so-called because of the apparent impossibility of escape from them 
- are not entirely black. That was the intriguing claim made by Hawking 00 
nearly forty years ago. Examining the behaviour of quantum fields in the vicinity of 
a black hole, he showed that, far from being emission-free, it should emit a steady 
flux of thermal radiation, with a temperature proportional to k, the gravitational 
field strength at the event horizon: 


k R T = 




he 3 


2irc 8 ttGM 


( 1 ) 


With this remarkable result, Hawking completed a thermodynamic treatment of 
black holes that had begun with Bekenstcin Ell] and has since been continued 
by many others, e.g., la ei 0 ei ei mg. This field brings together the normally 
disparate areas of gravity, quantum theory and thermodynamics; a glance at the 
various fundamental constants appearing in Eq. (|TJ) makes this fusion clear. Thirst 
for understanding of the underlying connections between these mighty realms of 
physics provides ample motivation for the study of what has come to be known as 
Hawking radiation. 

A pre-requisite for any such study must be the acknowledgement that Hawking 
radiation is not without its own problems, both practical and conceptual. On the 
practical side, the predicted temperature - at least in the gravitational context 
in which it was first derived - is virtually untestable. A solar mass black hole 
would, according to Eq. ([Tj) , have a temperature of about 10 _6 K - six orders 
of magnitude smaller than the temperature of the cosmic microwave background 
(CMB). Any radiation from the black hole would be drowned out by the CMB. 
Therefore, experimental verification of black hole radiation would seem to require an 
extremely light black hole, orders of magnitude lighter than the Sun. Such “micro” 
black holes may have formed early in the life of the universe m, or they may be 
formed today in the high-energy collisions at the Large Hadron Collider [21 FTB] - 
Such tiny black holes would, according to Eq. ([Tj) , have a very high temperature 
and quickly evaporate in a burst of radiation; however, if they do exist, they have 
so far escaped detection. 

Conceptually, there is the trans-Planckian problem [TEl IT5] - which has to do 
with the validity of the derivation of Hawking radiation. Let us briefly explain 
the problem here; it is discussed mathematically in §3.3 In most derivations, the 
spacetime is assumed to collapse, as in a star collapsing to form a black hole mm- 
Modes of the quantum vacuum are incident from infinity, propagating through the 
collapsing spacetime and out to infinity again, experiencing a gravitational redshift 
as they climb out of the ever-deepening gravitational well. The steady thermal flux 
seen at late times can be traced back to those vacuum modes which just manage to 
escape the event horizon, slowed and redshifted to greater and greater degrees. Thus, 
any low-frequency mode seen in the late-time thermal spectrum can be traced back 
to an incident vacuum mode of ever-increasing frequency; indeed, to an exponentially 
increasing frequency! The frequencies of these incident modes very quickly exceed 
the Planck scale [T5], widely believed to be a fundamental quantum limit. We cannot 
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justify the use of quantum field theory at such scales - yet it appears that Hawking 
radiation is dependent on the existence of these initial trans-Planckian frequencies in 
order to generate the final low frequencies at which the radiation should be observed 
rat Whatever is the correct physics, can we be sure that it will preserve Hawking 
radiation? 

These difficulties can be tackled by appealing to artificial event horizons, or 
physical systems possessing horizons analogous to those of gravitational black holes 
punuis]. This idea was first proposed by Unruh [[19] , who found that perturba¬ 
tions of a stationary background fluid flow behave just as a scalar field in Lorentzian 
spacetime [191 ED] • In particular, if the flow velocity crosses the speed of sound, the 
surface where it does so is entirely analogous to a black-hole event horizon, and on 
quantizing the perturbation field, one predicts analogue Hawking radiation in such a 
system. Of course, this model is subject to the same trans-Planckian problem as the 
gravitational one. In reality, however, the trans-Planckian problem is avoided by the 
ubiquitous phenomenon of dispersion pb|]. That is, the behaviour of waves changes 
at different scales by mechanisms which are better-understood than quantum grav¬ 
ity. For example, the discreteness of atoms or molecules places a fundamental limit 
on the wavelength of sound waves. It was later shown that, even after taking high- 
frequency dispersion into account, Hawking radiation is still predicted, with the same 
temperature as Hawking’s dispersionless model [] [21, [22]. The conclusion is that 
the trans-Planckian problem is a mathematical artifact, while Hawking radiation 
exists quite independently of the physics at the high-energy scale. This discovery 
has prompted a great deal of interest in a range of black-hole analogue systems: in 
Bose-Einstein condensates (BECs) (23] [Ml EH] [26] [27] , in ultracold fermions [28], in 
superfluid helium |29[ 3() [3TJ, in water [32], in electromagnetism [33] S3], and in 

These black- 


optics [35) 36] [37]. (See 112.2 for an overview of current research.) 
hole analogies might not teach us about quantum gravity directly, though they can 
demonstrate the ways in which the Hawking spectrum might change in response to 
new physics at the Planck scale. Perhaps more importantly, they offer a chance to 
study Hawking radiation as a general phenomenon related not so much to gravita¬ 
tion as to the restless nature of the quantum vacuum; any insight into the origin of 
this vacuum radiation [3HJ [39, 5Dj is to be welcomed. 

This tutorial aims to give an introduction to the theory of Hawking radiation in 
analogue systems. Given the wide variety of possible analogue systems, the tutorial 
aims to be as general as possible. However, analogue systems will obey different wave 
equations with a wide variety of dispersion relations. The framework is here laid 
mainly within the context of acoustic waves in fluids, in the spirit of Unruh’s original 
insight mm- This point of view has the advantage of simplicity, in the intuitive 
nature of the system, its general relativistic form (see [|2]) and its straightforward 
generalization to dispersive media (see [Q. 

The tutorial falls naturally into two distinct parts. Part [T] is concerned with the 
analogy between the black hole spacetime, as generally understood in the gravita¬ 
tional context, and the spacetime corresponding to a moving medium. The latter 


1 Assuming dispersion and the “surface gravity” are not too strong, as we shall see in later 
chapters. 
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point of view is more general, leading to the wealth of analogous systems briefly 
mentioned above. We then go on to derive Hawking radiation in this generalized 
spacetime, paying close attention to the important steps and ingredients in the 
derivation. In Part |TT] we generalize the theory further by allowing the medium to 
be dispersive, and so of greater applicability to experimentally realizable systems. 
Hawking radiation is re-derived, but only important differences from the previous 
derivation are highlighted. Since there is no exact analytic expression for the Hawk¬ 
ing spectrum in the presence of dispersion, we consider numerical and analytical 
techniques for its calculation. Finally, we examine the results of these calculations 
for a simple case. 
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Part I 

Dispersionless model 


2 The black hole spacetime as a moving medium 

2.1 The Schwarzschild metric 


General Relativity identifies gravity with curvature of spacetime, described by the 
spacetime metric ds 2 [TT]J • It has a unique spherically symmetric vacuum metric: 
the Schwarzschild metric, which in coordinates (ts,r, 0,(p) takes the form 

ds 2 = ^1 — —'j c 2 dt 2 s — ^1 — — j dr 2 — r 2 dfl 2 , (2) 


where dQ 2 = d6 2 + sin 2 6 d(p 2 is the angular line element, ts is the Schwarzschild 
time coordinate and rs = 2 GM/c 2 is the Schwarzschild radius. This describes, 
for example, the spacetime exterior to a star or planet with relatively low rotation 
rate. As r^/r —^ 0, the Schwarzschild metric approaches the flat Minkowski metric, 
so the coordinates (fs,r, 0,0) correspond to the usual spherical coordinates of flat 
spacetime for an observer at infinity. However, Eq. ([2]) contains two singularities, 
at r = 0 and r — rs- Since the Schwarzschild metric is valid only in vacuum, these 
singularities are relevant only when the entirety of the mass is confined to a radius 
smaller than rs, in which case it will inevitably collapse to a single point of infinite 
density at r = 0. Such objects are called black holes. The point r = 0 is a genuine 
singularity of Schwarzschild spacetime BH. and we shall not be concerned with it 
here. It is the surface r = rs - the event horizon - that is of interest to us. 

Let us briefly review the effects of the event horizon by examining light trajec¬ 
tories, or null curves, with ds 2 = 0. For simplicity, we shall consider only radial 
trajectories, so we also set dVt 2 = 0. This leaves us with a differential equation for 
the radial null curves: 



(3) 


Far from the Schwarzschild radius, where r rg, \dts/dr\ —» 1/c, so that light 
behaves just as it does in flat spacetime. However, as we approach the Schwarzschild 
radius, \dts/dr\ diverges in such a way that light takes longer and longer to travel 
any distance, and, if travelling towards the event horizon, can never reach it in a 
finite time ts- As far as the Schwarzschild time ts is concerned, the event horizon 
is infinitely far away. 


2.2 The Painleve-Gullstrand-Lemaitre metric 

Despite this peculiar behaviour, the event horizon is not a genuine singularity of 
Schwarzschild spacetime BD; it simply appears as such in the coordinates ( ts,r ), 
in terms of which the two regions separated by the event horizon are infinitely far 
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apart. This problem can be addressed by following in the footsteps of Painleve [02], 
Gnllstrand [53] and Lemaitre pH], dehning a new time 


( = is + 2 ^ + ^ln 

c c 


\J r ! r s 


V r / r s + 1 


which, when substituted in Eq. ([2]), yields the transformed metric 


ds 2 = c 2 dt 2 — 



(4) 


(5) 


In the coordinates (t,r,6,<p), the metric clearly has no singularity at r = rg ■ At 
constant t, the spatial metric is precisely that of Minkowski space, and the distance 
to the event horizon is always finite. Notice that, while the definition of t, Eq. 0 , is 
applicable only when r > rg, the metric of Eq. (|5| is easily extendable to all values 
of r greater than zero. The new coordinate t has opened up a previously inaccessible 
region of the spacetime. Keeping t fixed while decreasing r, we see from Eq. Q 
that, as we approach the Schwarzschild radius, we must have tg — * oo to compensate 
for the divergence of the logarithm. Thus, with respect to our original coordinates, 
the transformation to the coordinates (t, r ) is accompanied by an extension of the 
spacetime into the infinite future. Our discussion of light trajectories anticipated 
this: since an ingoing light ray approaches the horizon at an infinitely slower rate 
(with respect to tg), crossing the horizon requires tg —» oo. 

Let us take a moment to interpret the metric ([5]). Again, we shall consider 
only radial trajectories, setting dVt 2 = 0. The key point to note is that, if dr/dt = 
—C\Jrgjr, the metric reduces to ds 2 = c 2 dt 2 . Since this condition clearly maximises 
ds 2 , these trajectories are geodesics, and t measures proper time along them. It is 
as though space consists of a (Galilean) fluid |05j 06], flowing inwards with velocity 
~CyJrg/r to converge on the point r = 0. The geodesics just defined are those which 
are stationary with respect to this fluid. They define a locally inertial frame which 
we shall call the co-moving frame, and in this frame - i.e., with respect to the fluid - 
the speed of light is c. At the Schwarzschild radius, the fluid flows inward with speed 
c; anything that falls beneath this radius, no matter its velocity with respect to the 
spacetime fluid, will inexorably be dragged towards the centre at r = 0. This view 
is reinforced by looking at the radial null curves in the coordinates (t,r). Setting 
ds 2 = 0 and dVt 2 = 0, we find two possible trajectories for light: 



The forms of these trajectories near the point r = rg are shown in Figure [l] The 
second possibility represents rays propagating with the fluid; the total velocity is 
—c — cy/rg/r, the sum of the light’s velocity and the fluid’s velocity. It is perfectly 
regular at the horizon, and tilts over as it propagates, travelling faster and faster 
as it moves to ever smaller radii. The first possibility represents rays propagating 
against the fluid, having a total velocity of c — c^/Vgfr. This is not regular at 
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the horizon, nor should it be: there, the competing velocities of the fluid and the 
light exactly cancel, giving a total velocity of zero. Rays at higher radii will have a 
positive total velocity, and will eventually escape to infinity; rays at smaller radii will 
have a negative total velocity, unable to overcome the fluid flow, and will propagate 
inwards to r = 0. 

2.3 A general metric for moving media 

This analogy with a moving medium forms the basis of artificial black holes and 
event horizon analogies. We may simply replace —c^/rgjr in Eq. (jH) with the more 
general velocity profile V (x) to obtain, in 1 + 1-dimensional spacetime, the metric 

ds 2 = c 2 dt 2 — ( dx — V (x) dt ) 2 , (7) 

where now c is to be interpreted as the velocity with respect to the medium in 
question (not necessarily the speed of light). The co-moving frame, then, is the 
frame in which this medium is at rest, and in which the wave speed is exactly c. 
By contrast, the frame with coordinates x and t shall henceforth be called the lab 
frame. We shall always assume the medium to be left-moving in the lab frame, so 
that V < 0. 

From now on, we shall consider the general metric ([7]), which need not be grav¬ 
itational in origin. It may, for example, be applied to a system so far removed 
from astrophysics as a river flowing towards a waterfall M. so that the flow speed 
increases in the direction of flow, as illustrated in Figure [2} Imagine this river is 
populated by fish who can swim only up to a maximum speed c with respect to the 
water. Then the above metric suffices to describe the trajectories of fish in this river. 
Fish who are far from the waterfall, where the current is low, are free to swim around 
as they please, experiencing no significant resistance in either direction. However, as 
the current increases, there may be a point at which \V\ = c. As the fish approach 
this point, they will find it increasingly difficult to swim back upstream; passing 
this point, motion upstream is impossible, for the current is so strong that the fish, 
no matter how hard or in which direction they swim with respect to the water, 
are doomed to be swept over the waterfall. The point where |Vj = c is the event 
horizon, and the trajectories of fish swimming at exactly the speed c are analogous 
to the trajectories of light near a black hole horizon [^] 

2.4 White holes 

Now consider a slightly different scenario, portrayed in Figure [3j Here, water is 
flowing from a waterfall, so that its flow is initially very fast and slows as it travels. 
There is also a point here at which \V\ = c, but it behaves in a qualitatively different 
way from that in Figure [2] Fish far from the waterfall, where the current is low, 

2 Slightly more realistically, it is the trajectories of waves in the water - sound waves [1911201145] 
or even surface waves [III HU HU - which respect the black hole analogy. It is most common in 
the literature to make reference to sound waves, and for this reason the flow velocity is usually 
classified as subsonic (|V| < c) or supersonic (|V| > c). These terms are adopted in this tutorial. 
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may come and go as they please; but as they travel upstream, they will find it 
increasingly difficult to continue, and must come to a complete standstill exactly 
at the horizon. Whereas before the fish could not escape from the region beyond 
the horizon, now they find they cannot enter it. This is a white hole spacetime. In 
more precise terms, the supersonic flow (see footnote [2| leads away from the horizon 
in a black hole spacetime, but towards it in a white hole spacetime. Equivalently, 
we may think of the white hole as the time-reversed black hole: in the spacetime 
metric ([7]), we substitute V —» — V to transform between them. This one-to-one 
correspondence between black and white holes persists in the presence of dispersion 
(see, e.g., Appendix D of [50]). 

The white hole is of limited use in astrophysics. It is encoded in Schwarzschild 
spacetime, from which it is derived by taking the negative square roots in Eqs. 
0-(§- but it extends the spacetime into the infinite past. Therefore, its validity 
requires that the Schwarzschild metric is valid for t$ —> — oo, which is seldom the 
case, as black holes are believed to form from gravitational collapse. However, this 
is specific to the Schwarzschild case, and there is no fundamental restriction on the 
existence of white holes in the more general case of moving media. Indeed, many 
experimental setups achieve sub- or supersonic flow only over a limited region, so 
that a black hole-white hole pair is formed. White holes have thus found relevance 
in the field of analogue Hawking radiation. 

One point of particular interest and controversy is the stability or otherwise of 
spacetimes containing white holes. While black holes are always found to be stable 
PH, S2], white holes were found in [51] to be intrinsically unstable, and, conversely, 
to be stable in [53]. A systematic numerical study in [52] found that the stability 
or otherwise of white holes depends crucially on the boundary conditions imposed, 
though the nature of the physically appropriate boundary conditions remains for 
now a moot point (see also §5.2.4 of [IS] for a concise description of this problem 
and the issues involved). Another spacetime showing signs of instability is the black 
hole-white hole pair. Periodic geometries including such a pair are found to induce 
narrow instability “fingers’ 5 in the parameter space [231 EH El]. Better understood, 
however, is the dynamical instability of the black hole-white hole pair in a non¬ 
periodic configuration, in which the inner region acts as a resonant cavity: the 
so-called black hole laser [551 1551 [57, 55] . 

For simplicity, in the following we shall restrict our attention to black hole con¬ 
figurations, but white holes can be treated in an analogous fashion (see [59] for an 
analysis and comparison of black holes and white holes). 

2.5 Dispersion 

There is one caveat to the generalized metric ([7]) for moving media, but it is a serious 
one. Assuming for simplicity a massless field, then c, a constant parameter, is the 
speed of waves with respect to the medium no matter their frequency or wavevector. 
The medium thus described is dispersionless: the wave speed is absolutely fixed. In 
some ways, this is a good thing: the wave equation is simpler; the general solution is 
more tractable; the positions of any horizons are defined absolutely. But, as we shall 


soon see, it is precisely the absence of dispersion^ that leads to the trans-Planckian 
problem]^} So, while the generalized metric ([7]) may free us from the specific case of 
a gravitational black hole, it does not rid us of this fundamental conceptual issue. 
In the case of gravity, quantum theory is thought likely to modify things near the 
Planck scale, introducing high-frequency dispersive effects into General Relativity; 
but a quantum theory of gravity is as yet unknown. In real media that may be 
used in experiments, however, dispersive behaviour is generally well-understood, so 
that it may be possible to incorporate it into the theory of Hawking radiation. We 
examine how this can be done in Part El 


3 The wave equation and its solutions 


3.1 Deriving the wave equation 

Having now generalized the metric of a Schwarzschild black hole to that of a moving 
medium, let us study the behaviour of fields in a (1 + l)-dimensional spacetime 
described by the metric ([7]). For simplicity, we shall assume a massless scalar held; 
in the context of a real fluid, for example, such a held arises from small perturbations 
in the background how (see p3] for a detailed derivation). From the outset, we shall 
treat the held as complex, for even though we later (in ( |4.3 ) constrain it to be real, 
its decomposition into complex-valued modes is fundamental to the quantization of 
the held and, ultimately, to the Hawking process itself. 

Let us begin with the Principle of Least Action: the held (ft (t, x ) varies from one 
configuration to another in such a way that the action is an extremum (usually a 
minimum). The action is the integral 


5 = 



dx dt L (eft, (ft*, d t (ft, d t (ft*, d x (ft, d x (ft*), 


( 8 ) 


so the physics of the model is completely contained in the Lagrangian density, L. 
An extremum of the action is found by infinitesimally varying the helds (ft and (ft* 
and their derivatives, then setting the resulting variation in S to zero. This yields 
the Euler-Lagrange equation, 


dL d 


dL 


dt \d ( d t (ft* 


d 


dL 


dx \d (d x (ft*) 


= 0 . 


(9) 


The Lagrangian density for a massless scalar held is 


L = , 


( 10 ) 


3 More precisely, it is absence of high-frequency dispersion that is problematic. Thus massive 
fields, though they are dispersive at low frequencies and would allow us to retain the equivalence 
between medium and Lorentzian spacetime, do not resolve the trans-Planckian issue. In any 
case, they would restrict us to media with a very specific dispersion, whereas realistic media yield 
complicated dispersion relations and inevitably break the medium-spacetime equivalence. 

4 Although dispersive profiles are also problematic if the asymptotic velocity vanishes; see “Con¬ 


ceptual issues” off 12.2 and §VI of [B01 . 
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which is simply the covariant form of the corresponding Lagrangian density in flat 
space; is the metric tensor, g is its determinant, and g^, with raised indices, is 
its inverse. From the metric (J7|), we find that the Lagrangian density is 

L = ^{\{d t + v d x )(j)\ 2 - c 2 \d x (p\ 2 ) , (11) 

and plugging this into the Euler-Lagrange equation ([9]) yields the wave equation 

(d t + d x V) (d t + Vd x ) 0 - c 2 dl4) = 0 . (12) 

The partial derivatives act on everything to their right, including factors of V. 


3.2 General solution: the u- and ^-branches 


A general solution to the wave equation (12) is easily found. We define new variables 
u and v as follows: 


u = t — 


dx' 


c + V ( x ') ’ 


V — t + 


dx' 


c — V (x') 


(13) 


On substitution in the metric ([7]), we find that the metric in the coordinates (u,v) 
takes the form 

ds 2 = (c 2 — V 2 (x)) du dv , (14) 

which in turn, via the Lagrangian (10) and the Euler-Lagrange equation ([9]), leads 
to the wave equation 

d v d v <f> = 0. (15) 

So, in the absence of horizons, 0 is simply a sum of two arbitrary functions, one a 
function of u only, the other a function of v only: 


0 = 4>U (u) + 0u (v) 


t — 


dx' 


c + V (x') 


+ 


t + 


dx' 


c — V (x') 


(16) 


In the co-moving frame, (j) u is right-moving (counter-propagating - recall we take 
V < 0) while 0„ is left-moving (co-propagating). The fact that the two functional 
forms maintain their shapes is a consequence of the absence of dispersion: all wave 
components have the same velocity, c, with respect to the fluid, and so does the 
waveform as a whole. The only ambiguity is in the direction of travel; thus the 
solution splits into a right-moving and a left-moving part. 


Note that the metric (14) is related to the metric of a stationary medium (in 
which V is identically zero) simply by multiplication by a coordinate-dependent 
prefactor. This is a conformal transformation , and its usefulness stems from the 


invariance of the Lagrangian (10) - and consequently of the wave equation (15) 


under such a transformation. This observation corroborates a well-known theorem 
that all curved two-dimensional spaces are conformally fiat (i.e., related to flat space 
via a conformal transformation); see, e.g., 
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3.3 Event horizon 

In the vicinity of a horizon where V = —c, inspection of Eqs. ( Jl3| ) shows that u 
diverges as the horizon is approached. So the range u G (— 00 , 00 ) applies only to one 
side of the horizon, and the coordinate pair (it, v) does not cover the entire spacetime. 
To include the region on the other side of the horizon requires the introduction of an 
additional w-coordinate, and, consequently, an extra function in the general solution 

Let us suppose that the horizon is situated at the origin, x = 0; and, moreover, 
that the derivative of V is non-zero there. Then we may approximate the flow 
velocity to first-order in x as follows: 


V (x) ~ —c + ax. 


(17) 


where a > 0 for a black hole horizon. (A white hole horizon would have a < 0. 


Equation (16) states that counter-propagating waves are described by the arbi¬ 


trary function (j) u (u), where u is given by 


u = t 


dx' 


c+V (V) 


1 

t - 

a 


log (( |z|) , 


(18) 


the second equality holding in the vicinity of the horizon. It is evident from Eq. 


(18) that a horizon is located at x — 0, for the space is divided into two separate 


regions, with x = 0 marking the boundary between them. To see this, imagine 
first that we have a wavepacket centred at a certain value of u, and that this is 
located, at a certain time, at a positive value of x\ such wavepackets are illustrated 
in Figure [lj If t increases, then log (a \x\ /c) /a. must increase by exactly the same 
amount, and so |x| increases and the wavepacket moves to the right. Similarly, if 
we trace the wavepacket back in time by decreasing t, log (ct |x| /c) /a must also 
decrease by the same amount, and x decreases; the wavepacket has come from the 
left. But the logarithm diverges to —00 at the origin; this means that, no matter 
how far back in time we look, log [a\x\ / c) / a can be decreased by a corresponding 
amount without x ever becoming negative. The wavepacket must have originated 
arbitrarily close to the event horizon, moving very slowly forwards at first, picking 
up speed the further it travels. The longer the horizon has existed, the closer to 
it the wavepacket must have been at the moment of formation; furthermore, the 
thinner it must have been, since the lines of constant u bunch together ever closer 
as they approach the horizon. As we trace it back into the asymptotic past, a given 
wavepacket at a given position must have arisen from a wavepacket arbitrarily close 
to the horizon, of arbitrary thinness and composed of arbitrarily short wavelengths. 
This is precisely the trans-Planckian problem. 

Exactly the same analysis holds in the region to the left of the horizon, as also 
shown in Fig. 06)- 

We shall always take V to be constant in time, so that, mathematically at least, 
any horizons must have existed since the infinite past (though this is a subtle issue, 


and we shall return to it in §15.1 and 6.2). In this case, all wavepackets, when traced 


backwards in time, are found to have originated arbitrarily close to the horizon and 
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can never have been in the region on the opposite side. Since a complete description 
of the field must apply throughout the entire space, we conclude that there must 
be two distinct w-coordinates, each localized to one side of the event horizon. Let 
us denote them Ur and Ul- In the definition (13), these are distinguished by the 
lower bounds of the integral over x, that for Ur being in the right-hand region and 
that for Ul in the left-hand region. Then the coordinates Ur and Ul both range 
from — oo to +oo, but cover only their respective halves of the spacetime, as seen 
in Fig. [lj We thus see that the existence of a horizon - or, equivalently here, of 
one subsonic and one supersonic region - results in three independent solutions of 
the wave equation ( Ur , Ul and v ) rather than the usual two (u and v). Ur and 
Ul are both of w-character since both describe waves travelling to the right in the 
co-moving frame of the fluid. The difference between them is in the lab frame: ul 
applying in the supersonic region where the fluid flow is faster than the wave speed, 
the ul- waves are seen to be dragged to the left in the lab frame; whereas, since ur 
applies in the subsonic region where the fluid flow is less than the wave speed, the 
Ur- waves still manage to move to the right in the lab frame, albeit at a reduced 
speed. 

The requirement of two w-coordinates in addition to a single u-coordinate to com¬ 
pletely describe the spacetime answers an immediate objection that might be raised 
against the possibility of Hawking radiation in 1 + 1-dimensional spacetime: namely, 
the conformal flatness of all two-dimensional spacetimes, mentioned previously in 
§3.2 At first sight, this seems to suggest that all two-dimensional spacetimes are 
conformally equivalent, thus precluding the occurrence of Hawking radiation in some 
of them (those with horizons) given its absence in others (without horizons). But 
conformal flatness does not imply conformal equivalence if there exist non-equivalent 
spacetimes which are flat nonetheless. This is precisely the case here: whereas the 
horizonless spacetime is conformally equivalent to the two-dimensional plane, the 
existence of the horizon splits the spacetime into two regions, each of which is 
conformally equivalent to the two-dimensional plane, these planes being “stitched 
together” at the horizon [61]. Therefore, the flat equivalents of spacetimes with and 
without horizons are topologically distinct; they are not conformally equivalent to 
each other, and the absence of Hawking radiation in one does not imply its absence 
in the other. 


4 Field modes and quantization 

4.1 Stationary modes and the dispersion relation 

The Lagrangian density ( [II] ) is invariant under time translation; that is, performing 
a small shift in the time coordinate, t —> t + At, leaves L unaffected. This implies 
[63 ] the existence of stationary modes of the form 

(f(t,x) = e~ lU}t (f) u {x) . (19) 


On quantization, these correspond to energy eigenstates or quasiparticles 


Equa¬ 


tion (16) indicates that, in the absence of an event horizon, the spatial part of the 
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solution takes the form 


(x) = C u exp ( iui 


dx' 


c + V ( x') 

= c v <%(x) + c v <i>z(x), 


+ C v exp ( — iu> 


dx' 


c — V (x') 


( 20 ) 


whereas, accounting for the considerations of §3.3[ in the presence of an event horizon 
(at x = xh) we have instead 


(x) = C UR exp ( iu 


dx' 


'XR 


c + V (x') 


9 (x — xh) + C UL exp ( iu 


dx' 


’X L 


c + V (x') 


9 (x H - x) 


+C V exp ( —iuj 


dx' 


c — V (. x') 


= C UR (x) 9{x- x H ) + C UL (pl L (x) 9 (. x H - x) 

+C V ft (x) , 

where ffeaviside step functions (9(x) = 0 for x < 0 and 1 for x > 0) have been 
included to clarify the restricted domains of and Like the general solution 
(16), the general stationary mode for a given frequency splits into two independent 
parts, one counter-propagating and one co-propagating [^J which we henceforth refer 
to as u- and w-modes, respectively; and, in the presence of a horizon, the w-branch 
splits further into modes localised to either the subsonic or supersonic region. 

Defining the local wavevector k (x) such that <0 V (x) = exp f x dx' kd v (x') j, 
we see that 


K (*) = 


UJ 


K {x) = 


LU 


c + V (x) ’ 

relations summarized by the single equation 

(u> — Vk) 2 = c 2 k 2 


V(x) ’ 


( 22 ) 


(23) 


( 21 ) 


This formula is very intuitive once we recognise that the lab frequency oo is related 
to the co-moving frequency (n cm via the Doppler formula: c u = (u cm + Vk^\ Equation 
(23), then, is simply oj 2 m = c 2 k 2 - the dispersion relation in the co-moving frame. 
Comparing with Eqs. (22), we see that 


W cm = UJ — Vk 


ck for u— modes 
—ck for v— modes 


(24) 


5 Being stationary modes, the wave envelope does not move at all. A way of thinking about the 
wave velocity is to consider a wavepacket, strongly peaked at the given values of u> and k. Then the 
wavepackets corresponding to the u- and v-modes will travel in opposite directions with respect to 
the fluid. Wavepackets are included in Fig. |T| 

6 The “spacetime fluid” itself is Galilean with respect to the lab frame - the spatial coordinate 
is transformed according to dx' = dx — Vdt in Eq. 0 - which is why the Doppler formula appears 
in Galiliean form. This is in contrast to relativistic motion occurring within the fluid [JB]. 
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The function ±cfc will be referred to as the dispersion profile or curve. In this 
dispersionless case it is simply proportional to k, but general dispersion profiles are 
more complicated, and will be considered in Part [TTJ 

Equation (24) - which states that the co-moving frequency must both satisfy the 
dispersion relation in the co-moving frame and be related to the lab frequency via 
the Doppler formula — can be solved graphically from one of two possible viewpoints: 


• Co-MOVING FRAME: The dispersion profile in the rest frame of the fluid is 
plotted alongside the co-moving frequency, which is given by the straight line 
U) — Vk. The points of intersection occur at the possible wavevector solutions. 
In this picture, variation of the velocity V does not alter the dispersion prohle, 
but the slope of the Doppler curve, whose ^-intercept is equal to uj and is 
therefore to be kept constant. An illustration, with V varying between sub- 
and supersonic values, is given in Figure |4j 

• Lab FRAME: The lab frequency ui, which is constant and appears as a straight 
horizontal line, is plotted alongside Vk ± ck, the dispersion prohle as viewed 
from the lab frame in which the fluid is not at rest. Again, the points of 
intersection give the possible wavevector solutions. In this picture, variation 
of V alters the dispersion curve, causing it to tilt on a dispersion diagram. 
The equivalent of Figure [4] as viewed from the lab frame is shown in Figure [5} 


4.2 Scalar product 

The Lagrangian L is also invariant under phase rotation: 0 —y 0e*“, where a is 


a real constant. For any two solutions and 02 of the wave equation (12), this 
symmetry in the phase implies [53] conservation of a scalar quantity - the scalar 
product (52] - defined as 


P+OO 


( 01 , 02 ) = i I dx { 0 * (d t + Vd x ) 02 02 (d t + Vd x ) 0 *} 

J —oo 

/ +oo 

dx{4>\Tl2 - 02 <} , 

-oo 

where the canonical momentum is 


dL 


IT = 


d (d t <j>*) 


= (dt + Vd x )<p. 


(25) 


(26) 


The definition (|25|) of the scalar product directly implies the following relation be- 

(27) 


tween complex conjugate solutions: 

(01,0a) = -(01,02)* • 


The scalar product of a solution with itself is called the norm. Note that this norm 


is not positive-definite: if 0 solves the wave equation (12) and has positive norm, 
then the complex-conjugate solution 0* - which also solves the real wave equation 


(12) - must, according to Eq. (27), have negative norm 
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The scalar product (25) only makes sense for a complex field, vanishing when 
0 is real. However, even when the total field is real, it makes sense to decompose 
it into complex components, and these will evolve in such a way that their scalar 
products are conserved. 


4.3 Field decomposition: the ^-representation 

Quantization of the field requires expressing it as a sum over eigenstates, which are 
then quantized individually. For now, let us assume that V is constant in x: then 
L is invariant under spatial translations and k is also a conserved quantity [63] . 
We may then work with the states of constant k, which correspond to momentum 
eigenstates, and for which Fourier analysis assures completeness on the entire real 
line. This might seem a trivial case, but it will turn out to be very instructive for 
the generalization to inhomogeneous velocity profiles. 


4.3.1 Orthonormal basis 


The canonical momentum operator (d t + Vd x ), defined in Eq. (26), is the derivative 
with respect to time in the co-moving frame: when acting on a stationary mode, 
it multiplies it by —iu> cm = ^fick, where the minus (plus) sign corresponds to 
u- (■ v -) modes. V being constant, the stationary modes are simply plane waves 


exp (ikx — iut), whose scalar products (see Eq. (25)) are 

s e ik 2 x-iu> u (k 2 )t^ = ^ c ^ g ^ ? 

z ik lX -i^{ kl )t , e ik 2 x-i^{k,)t^ = c k 1 S(k 1 ^k 2 ) , 


I 


0 ik\x—iu u (k\)t ik2X—iuj l 


( fe 2 )t\ _ Q 


(28) 


(Since V is constant, there is no event horizon and hence no splitting of the u- 
branch into Ur and Ur.) That is, the plane waves are mutually orthogonal, and, 
utilizing Eqs. (24), their norms are 47ra; cm 6 {k\ — fc 2 ). Although we can renormalize 


these modes using an appropriate prefactor, the occurrence of one mode and one 


complex-conjugate mode in the definition (25) of the scalar product ensures that we 


can never change the sign of the norm. Therefore, the sign of the norm is equal to 


the sign of the co-moving frequency, and the dispersion relation (23) leads also to 
the classification 


Wcm = UJ — Vk = 


\ck\ for positive — norm modes 
— \ck\ for negative — norm modes 


(29) 


Combining classifications (24) and (29), we see that positive-norm modes are those 
which: 


• if k > 0, are u-modes; 

• if k < 0, are rmnodes; 
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while for negative-norm modes, this correspondence is reversed. (In Figs, [f] and 
[5j the splitting of the dispersion profile into the u- and u-branches, and into the 
positive- and negative-norm branches, is labelled explicitly.) Thus positive-norm 
modes correspond to the intuitive notion that the sign of k determines the direction 
of travel (with respect to the medium); and the negative-norm modes, through 
complex conjugation, correspond to precisely the opposite. 

The orthonormal A;-mode basis, then, is the set of modes 

0^ {x, t) = 1 exp (ikx - iu u (k)t) 

A/47T \ck\ 

<t>l (■ x , t) = 1 exp (ikx - iuj v (k)t ) 

a/47t \ck\ 

normalized to ± S(k—k') where the sign is determined as above, and where uj UJ,v (k) = 
(V ± c)k. 


= 0“ (x) e~ iojU ( k)t , 

= 0l (x) e -‘^‘, (30) 


4.3.2 Field decomposition 

According to Fourier analysis, the set of all real wavevectors form a complete set in 
position-space. The separate u- and w-branches found here are needed to describe 


the full time-dependence of the field, the wave equation (12) being second-order in 


time. So the general solution of the wave equation can be written in the form 


*+00 


0 (x, t)= dk [a u (k)4% (x) e ~ iujU{k)t + a v (k)(f>l ( x ) e ? ( 31 ) 


where a u (k) and a v (k) are complex-valued functions of k. It is useful to separate 
the positive- and negative-norm modes in the integral, as follows: 

ro 

0 (x, t) = dk [a v (k)(pl (x) + a v *(k)(t>l* (x) e iujV{ - k)t \ 

J — OO 

poo 

+ dk [a u {k)(t) u k (x) e~ iu}U(k)t + a u *(k)(j) k * (x) (32) 

Jo 

where, for simplicity, the coefficients of complex-conjugate modes have themselves 
been written as complex conjugates so that the total field is real. The first term in 
each of the integrands, written without a *, is a positive-norm mode; its complex 
conjugate, written with a *, has negative norm. If the total field 0 is known, then 
the orthonormality properties of the modes allow extraction of the coefficients using 
the scalar product: 


a u (k) = (0jf, 0) , a u *(k) = — (0^*, 0) , 

«•(*) = (0*,0) , aV *(k) = -(0r,0) . (33) 


4.3.3 Quantization 


In the form (32), the field is readily quantized via the usual methods of Quan¬ 
tum Field Theory (QFT) [52] , We promote the real-valued field variable 0 to the 
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Hermitian operator 0, the coefficients a u ^ v {k) to operators a^ v and their complex 
conjugates [a u / v (k)]* to the Hermitian conjugate operators 


4>(x,t) = 


dk d v k (j) v k 


CLuQu [x)e 

roo 

+ dk 
Jo 




a k (p k (x) e - iu ^ k)t + (x) 


(34) 


Similarly, using the definition (26), the canonical momentum operator is 
r o 


TT (x,t) = 


dk 


d v k 7T v k (x) e - iujV{k)t + af 7T v k * (x) e iujV{ ~ k)t 

a u k Ttl (x) e~ iuU{ - k)t + af (x) e iujU{k)t 


/»oo 

+ / dk 
Jo 


(35) 


The relations (33) may be similarly quantized, yielding expressions for the mode 


operators in terms of the full field operators: 


f+OO 


a h = i 




a k ~ 


dx{r k ‘ M OM) - <* , 

30 r . -s'. 

dx 1 4> k (x) e~ lu>u< ^ t; k (x, t ) — 7T k (x) e~ luU ^ t (j) (x, t)| , (36) 


with analogous relations for the x-mode operators. The quantization procedure is 
completed by imposing the canonical commutation relations Q 


0 (X, t) , 7T (x', t) =l5(x — x'), 0 (X, t) , 0 (x', t) = 7T (x, t) , 7T (x', t) 


= 0. 
(37) 


Finally, from these it can be shown that the mode operators satisfy the Bose com¬ 
mutation relations 

= 5 (k - k') , (38) 

all other commutators being zero. Thus we see that a k v and - the quan¬ 
tum amplitude operators multiplying modes of positive and negative norm - are 
bosonic annihilation and creation operators, respectively, for the u- or rmnode with 
wavevector k. 

This last point is of paramount importance, and worth emphasizing: positive- 
and negative-norm modes correspond to bosonic annihilation and creation opera¬ 
tors, respectivel}|^} Notice that, although the notation was chosen with foresight, 
we did not make any mathematical assumptions about the mode operators. Their 

7 We work in natural units where % = 1. 

8 Had we considered a complex rather than a real field, the field operator 0 would not be 
Hermitian and the annihilation and creation operators appearing in its decomposition would not 
be Hermitian conjugates of each other. They would then correspond to different particles, i.e., 
particles and antiparticles. Here, with a Hermitian field operator, we might say that we are 
considering particles which are their own antiparticles. 


n u 


n v 

a k , a k , 


a k , a k . 
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bosonic property (38) follows from the norms of the modes they multiply and the im¬ 
position of the canonical commutation relations (37). The correspondence between 
norm and bosonic quality represents a generalization to moving media of QFT in 
flat space. There, it is standard to separate the complex-conjugate positive- and 
negative-frequency modes, finding as here that the corresponding mode operators 
obey the Bose commutation relations (38). What we have just shown is that it is 
the sign of the norm which determines these mode operators, and this happens to 
be the same as the sign of the frequency in the co-moving frame - not the sign of 
the frequency in the lab frame. Worded this way, this point seems obvious, since for 
the constant-velocity profile we could have transformed into the co-moving frame, 
in which case the co-moving and lab frequencies would have coincided and flat-space 
QFT would have been applicable. Ultimately, however, it is the norm of the modes 
that enters the derivation of Eqs. (38), and this will continue to be true for inho¬ 
mogeneous profiles, where the co-moving frequency is no longer so well-defined and 
the norm emerges as the true fundamental quantity. 

The quantization procedure above is admittedly rather abstract. An analogous 
procedure, from a physically grounded viewpoint, is performed in condensed matter 
theory using the Bogoliubov theory of dilute Bose gases; see, for example, reference 


4.4 Field decomposition: the (^-representation 


Plane waves (30) of constant k are valid when V is constant - even if it is not 


constant everywhere. They are especially useful in the const ant-velocity asymptotic 
regions, where they correspond to quasiparticles that can be measured at infinity. 
An exact solution for a stationary inhomogeneous velocity profile will be expressible 
as a sum over plane waves - of equal lab frequency ui - in the asymptotic regions. 
This can be interpreted as a mixing or coupling between the various plane waves. 

Therefore, it is much more useful to consider oj as the integration variable when 
summing the modes [59], since this groups together all plane waves that can be mixed 
together by an inhomogeneous flow. First, we transform from the /^-representation 
to the cu-representation when V is constant; we then consider the generalization to 
inhomogeneous velocity profiles. 


4.4.1 Homogeneous flow 

For reasons that will become clear, let us begin with the ^-representation that does 
not separate positive- and negative-norm modes: 


i(x,t) = 


(39) 


where, consistently with Eqs. (31) and (32), we have (x) = 4>f* (x) and af k = 
and similar identities for the u-modes. Recall from £4.3 that the u -modes with 


k J 


k > 0 and the u-modes with k < 0 have positive norm. 
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We wish to express the field operator (39) in the equivalent form 


i (x, t) 



duj [al4,l(x)e- <ut 


+ (x) e-”‘] 


( 40 ) 


This requires a redefining of the modes and their corresponding operators. Firstly, 
we note that, in order for the quantization procedure to carry over exactly, the 
^-representation of the modes and operators should be normalized with respect to 

id: 


(C>C0 = S -w') 




= 6 (id 


id 


(41) 


and similarly for the w-modes. Now, the Dirac 5 function 6 (id — id') = \dk/du \ 5 (k — k'), 
while the differential did = \did/dk\ dk. So, if the modes and operators in the k- and 
(^-representations are related via [59] 


dk 1 


did 

dk u 


4 ) k u (uj) i 0OJ 


dk v 


did 


(pk v (uj) 1 


a u = 


did 


dk u (ui) 


dk 1 


a,., = 


did 


dk v (uj) 


(42) 


(43) 


we find that the normalization conditions (41) are satisfied and that the integrand 


of the u -representation of the field (40) transforms exactly into the integrand of the 


^-representation (39), since the factors of \did/dk\ cancel. So the orthonormal basis 


in the (^-representation is the set of modes 

C(M) = 1 

4>Z ( x , t) = 


—iujt 


\J 477 | ck u (id) Vg ( k u (ld )) 

1 

\J 47T I ck v (id) Vg ( k v (ld )) 


exp (ik u (id)x — iidt) = 0“ (x) e 


exp (ik v (id)x - iidt) = C (®) e , (44) 


normalized to ± 6 (id — id') according to the sign of the norm. We have defined 
v g (k) = dd/dk , which equals uj/k in the absence of dispersion, so that k(id) v g (k(id)) 
appearing in Eqs. ( [44] ) is simply equal to id. (In the dispersive case, v g is the group 
velocity; see ©) 


As before, we would like to separate the (n-modes of the field operator (40) into 


positive- and negative-norm modes, i.e., into those terms corresponding to annihi¬ 
lation and creation of quasiparticles. We have already seen in §4.3| that the sign of 
the norm is simply related to the wavevector k - but there is a non-trivial aspect 
in its relation to frequency. This can be clearly seen in the dispersion diagrams of 
Figs. [4] and [5] For d > 0. the solution on the u-branch always has positive norm; 
the u-rnode, on the other hand, has positive norm when the flow is subsonic, and 
negative norm when the flow is supersonic. Restricting our attention to the u-part 
of the field, we have 


V (x, t) = 


f S^did 

(x) e lU}t + a“V“* ( x ) e lut 


[/o 00 ^ 

(■ x) e iut + a u \(t) u * u (:c) e 

-iut 


for \V\ < c 
for \V\ > c 


(45) 
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where, as before, the first terms of the integrands are the positive-norm modes 
and their complex conjugates have negative norm. We emphasize the crucial fact 
that, comparing subsonic with supersonic flow, w-modes of equal lab frequency have 
opposite norm. 


4.4.2 Inhomogeneous flow 

There is nothing particularly significant about the sign of the frequency when the 
flow is truly homogeneous, for we could easily shift between subsonic and supersonic 
flow via a simple change of reference frame. There is no gain in transforming to the 
^-representation, for k itself is a conserved quantity, and mode mixing cannot occur. 

For an inhomogeneous flow, however, the significance is profound. There is a 
preferred reference frame in which the flow profile is time-independent, and in which 
oj is conserved. It is in this frame that the ^-representation is useful, in which the 
linking of modes with equal frequency has physical meaning. And if, in this frame, 
the asymptotic flow velocities are one subsonic and the other supersonic - if, in other 
words, an event horizon is present — then, through the switching of norm induced by 
this transition, positive-norm modes on one side can couple to negative-norm modes 
on the other. 

We know the form of the u-modes from 94.11 and we have seen in 93.31 that we 


must include two distinct modes, one corresponding to each side of the event horizon. 
Taking account of the norm switching as V varies between sub- and supersonic 
values, the spatially distinct positive-norm w-modes can, for uj > 0, be written 


0, t) = 9 0) - exp (-iuu R ) , 

V47rca; 

= 6 (~ x ) -J== exp (iuu L ) , 


a/4 


7 TCU 


(46) 


where Ur and Ul are the w-variables defined on the right- and left-hand sides (i.e., 
the subsonic and supersonic regions), respectively. Since those of the same frequency 


can mix, we can form linear combinations of one of the modes in Eqs. (46) with the 
complex conjugate of the other. To normalize such a linear combination, we simply 
employ the linearity of the scalar product: 


,L> a 4k, 2) R + /30-^ 2 ,Z,) 

= (c 1)R , 0“ 2 , r )+w/3 (c 1)R , ci 2 ,l) +r« (<r wl ,L, 


(M 2 - |^| 2 ) 5 (uq - uj 2 ) ■ 


So, if |ck | 2 — \/3\ 2 = 1, the linear combination R + L is normalized with 
positive norm; it is also automatically orthogonal to any such combination involving 
the modes L and 4^ri or modes of a different frequency, since the individual 
modes are orthogonal amongst themselves. Therefore, any set of modes 



ol u,\4 > 'L,r 


+ @W, 10 


U-k 

—lj,L ? 



«u, 20 


u 

—u),L 


+ Pu, 20 


U-k 

LJ,R 1 


(47) 
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where 


K/-|AJ 2 = 1, (48) 

form a complete set of orthonormal, positive-norm ?j- modes. The various possibilities 
correspond to different ways of characterizing the quasiparticles, and they are made 
possible by the horizon’s splitting of the w-branch into two independent parts. 


5 Hawking radiation 

5.1 In- and out-modes 


Despite the limitless possibilities that the freedom in the choice of u-modes allows, 
very few of them are useful. We should deal only with “natural” sets of modes that 
correspond to possible measurements. Consider the simplest example: a u j = 1, 
— 0, which is simply the set of localized modes given in Eqs. (46). These 
correspond to single outgoing wavepackets. A quasiparticle travelling to the right 
with frequency u is precisely an excitation of the mode 0“ R ] similarly, a quasiparticle 
travelling to the left with frequency —oj is an excitation of 0“ wi . They differ in 
this respect from all other linear combinations, which correspond to two outgoing 
wavepackets rather than just one. In light of this property, the set of modes 0“ R 
and 0“ wi are termed out-modes : those which correspond to a single outgoing wave 
in the asymptotic future. The w-modes 0" also have this property, and are therefore 
the out-modes of the w-branch. 

Complementing the out-modes are the in-modes , those corresponding to a single 
ingoing wave in the asymptotic past. The v -modes 0" are such modes, so that, 
on the u-branch, the in-modes are equal to the out-modes; this is a consequence of 
the regular behaviour of the u-coordinate at the horizon. In contrast, as we saw in 
1 3.3 the w-modes, in both the subsonic and supersonic regions, drift away from the 
horizon. How, then, can we possibly form an ingoing w-mode? One possible answer 


(see 16.2 for another) is that, in the asymptotic past, we are able to form ingoing 
waves if the event horizon is not present there. That is, we may assume that the 
horizon has not been forever present, but was formed at some instant in the past, 
before which the flow was everywhere sub- or supersonic. In this initial, horizonless 
spacetime, the u-modes are well-behaved and are precisely the in-modes we are 
looking for |^j There being no horizon at this initial stage, there is no divergence 
of the wavevector at any point, so that - although it will vary from point to point 
if V is inhomogeneous - the sign of the wavevector can never change. Positive- 
norm in-modes are thus composed entirely of positive wavevectors, and similarly 
the negative-norm modes of negative wavevectors. 

We are here interested only in the long-time, stationary state brought about by 
the final presence of the horizon, and not in any transient effects reliant on the details 
of its formation. It is convenient to refer only to the final spacetime with horizon 
present, and so, to find the form of the in-modes that must be used, they must first be 


°The it-modes of the horizonless spacetime would also be out-modes if this spacetime were not 
to form an event horizon; but, since the out-modes are by nature defined with respect to the 
asymptotic future, the true out-modes are those corresponding to the presence of the horizon. 
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propagated from the initial spacetime to the final one. Since the precise evolution is 
unimportant, we may consider simply an instantaneous change of the velocity profile. 
Immediately after this change, the spatial form of the in-modes is precisely as in 
the initial horizonless spacetime - but they will no longer be stationary. We may, 
however, form linear combinations of them to find stationary in-modes appropriate 
to the final spacetime, so long as we form such combinations only from modes 
of the same norm ; the reason for this is made clear in 15.2 As already noted, 


this requires that positive-norm modes be formed only from positive wavevectors. 
This implies in turn that the positive-norm modes, if analytically continued onto 
the complex x-plane, must be analytic in the upper half of this plane mmm- 


Recalling from Eq. (18) that, in the vicinity of the horizon, the coordinates Ur and 


Ul go as log ((a/c)|x|) /a, traversing the horizon x — 0 on the upper-half x-plane 
analytically connects these two coordinates with the addition of an imaginary part: 


Ur —y Ul -\- i'K /ol and Ur —>■ Ur — in/a. Applying these continuations to Eqs. (46), 
we find that the u-modes on either side of the horizon also analytically connect 
to each other, but, upon exponentiation, the relative imaginary component of the 
coordinates translates into a relative amplitude between the right- and left-hand u- 
modes. This relative amplitude determines the a and (3 coefficients of Eqs. (47), and 


after normalizing according to Eq. (48), we fold that the stationary positive-norm 


in-modes of the w-branch are, for ui > 0, 


iu, in 

Vr 


,u,m 


A sinh W) 

1 

y / 2 sinh (w) 


I 5^ ,u,out . — 5^ ill, out* \ 

\e 2a( t>J,R +e 2 “0-b,L ) > 

l 5^ iU.OUt I — / u,out*\ 

( e “"P4.,z, + e ) > 


(49) 


where f\f°R and are simply the localised modes of Eqs. (46) with their out- 

mode character made explicit. The negative-norm in-modes are simply the complex 
conjugates <f>j )R and 0J WiL . 

Despite the mathematical elegance of this derivation of the stationary in-modes, 
it is physically obscure. It is therefore of value to consider briefly the evolution of 
wavepackets, as is usually done in more direct physical treatments of the problem 
HI El US] - the stationary modes of Eqs. (091) may then be thought of as linear 
combinations of wavepackets. Those wavepackets located far from the horizon at 
the time of its formation are essentially unaffected by its appearance, except for a 
possible change of direction due to the transition between subsonic and supersonic 
flow. This accounts for the regularity of the in-modes away from the horizon. For 
those wavepackets located near the horizon at the time of its formation, there will be 
a marked shift in velocity and a marked change in its wavevector as it slowly moves 
away from the horizon, so that the in-modes exhibit the divergence in wavevector at 
the horizon characteristic of the n-modes. The wavepackets we have so far considered 
evolve into outgoing wavepackets on one side of the horizon, and therefore behave 
just as the out-modes do. The crucial difference, then, for the in-modes is that there 
are wavepackets which actually cross the horizon at the time of its formation. Such 
wavepackets will then evolve as two “disconnected” pieces - but they are not quite 
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disconnected, because in the past they are analytically connected to each other. 
This is the origin of the analytic continuation between the right- and left-hand sides 
across the horizon, which is characteristic of the in-modes. Given the logarithmic 
divergence of the phase at the horizon, this analytic continuation can only be of the 


form (49) or their complex conjugates, depending on whether analyticity is imposed 


on the upper or lower half-plane. Straightforward algebra then shows that these two 
possibilities correspond to the positive- and negative-norm modes, respectively. 

In-modes and out-modes form two distinct and natural sets of orthonormal u- 
modes, and we are free to decompose the w-part of the held operator in terms of 
either. Together with the t’-modes, they form a complete set of orthonormal modes 
that solve the wave equation (12). Therefore, the total held operator may be written: 




duj {Cy?C,!R (*, x ) + (*, x ) + d l4>l (*, x ) + h - c - 

duj {“uTffiT (*’ (*> x ) + d X (*> x ) + h - c -} ( 5 °) 


where h.c. stands for Hermitian conjugate, containing the negative-norm modes and 


creation operators. Substituting the transformations (49) into the expressions (50) 


for the total held operator, we find the corresponding transformation between the 
mode operators: 


-■ u, out 
l w,R 


OUt 

l -u,L 


sfi sinh (^) 


2 sinh (^) 


~«,in i ^u.inf 

e2 “ a J,R + e 2a a -u,,LI > 


( ~u,in , -5“ -w,int\ 

[ e 2 aa -u,L + e 2 aa J,R ) 


(51) 


It stands to reason that, since the modes of the in- or out-bases are formed from 
both positive- and negative-norm modes of the other, so the mode operators of one 
basis should combine both annihilation and creation operators of the other. 


5.2 Spontaneous creation 

Inequivalence of incoming and outgoing modes is not an exotic phenomenon in 
and of itself. Any scattering process exhibits this kind of behaviour: an incoming 
wavepacket or particle is dehected, often into several other wavepackets of particles, 
each of which can be considered a single outgoing mode. Normally, however, this is a 
process of conversion-, the incoming particles are converted into outgoing particles, 
such that the sum of the rates of outgoing particles - all contributing with the same 
sign - is equal to the rate of ingoing particles. Decrease the rate of incoming particles 
to zero, and the rate of outgoing particles will likewise vanish. We might say that 
all these particles have the same sign of norm. 

The difference here is that the converted particles have different signs of norm, 
leading to the normalization condition \a\ 2 — |/3| J = 1. Since the norm is always 
conserved, a part-conversion of the incoming modes into negative-norm modes must 
result in an increase in the amount of positive-norm modes. This process is, at least 
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in part, an amplification [HSEBj. For oppositely-normed particles equal numbers of 
both must be added to the system to keep the overall norm constant. Moreover, at 
the quantum level, an amplifier adds particles to the system even when there are no 
incident particles m- 

Mathematically, the vacuum state is that in which no modes are excited: it is the 
absence of quasiparticles. Since one cannot annihilate an excitation from the vacuum 
state, it must vanish when acted on by any annihilation operator. The vacuum is 
thus defined as the eigenstate of all annihilation operators with eigenvalue zero [62j: 

a |0) = 0 Vo. (52) 


Clearly, this depends on the particular set of annihilation operators, and hence on 
the particular basis of modes [68]; in particular, we can define both an in-vacuum 
and an oid-vacuum, according to whether the vacuum state has zero eigenvalue for 
all in-annihilation operators of out-annihilation operators. That said, if a mode 
transformation only mixes modes of the same norm, then the corresponding trans¬ 
formation between the mode operators will expand annihilation operators only in 
terms of other annihilation operators, and the state with zero eigenvalue for one 
set of these operators will also have zero eigenvalue with respect to the other: their 
vacuum states are identical. This is why, when forming the stationary final-state 
in-modes in (5.1, we could only combine modes of equal norm: this ensures that 
the in-vacuum state defined by the late-time in-modes is exactly the same as that 
defined by the initial in-modes. Similarly, if positive-norm in-modes were found to 
scatter only into positive-norm out-modes, then the vacua defined by the in- and 
out-modes would be identical, and an absence of incoming particles would lead to 
an absence of outgoing particles. 

On the other hand, if scattering into opposite-norm modes, and hence mixing 
of annihilation and creation operators, takes place - as it clearly does in Eqs. (49) 
and (51) - then the in-vacuum is not equal to the out-vacuum, and the absence of 
incoming particles must lead to some presence of outgoing particles [68j! This is the 


mysterious effect of spontaneous creation, which is purely quantum mechanical in 
origin 


10 


5.3 Radiation from an event horizon 


We may demonstrate this particle creation explicitly using Eq. (51). Imposing the 
in-vacuum, the expectation value of right-moving outgoing quasiparticles is found 


10 An analogous kind of creation takes place in the Unruh effect where the vacuum seen by 
an inertial observer is found to differ from the vacuum seen by an accelerating observer, leading 
to the detection of a thermal bath by the latter. 


24 








to be 


/ l£,OUt\ 

( n V ) = 


(o 

(o 


in I a uj,R 


U,OUtf 


lu'R 


|0,»> 


\/ 2 Sinh('“) 


^ ~u, int , 

e2 “Vi? +e 2a V 


x 




sy. ^u in . _2y-~u,int 

e 2 “V,fl + e 2 “V',l 


|0in) 


y/4 sinh (^) sinh (^) 

1 

7 e 

y/4 sinh sinh (^) 

1 -™ c / /\ 

-t— r e “ o leu — u ) 

2 sinh 

“W—- 6(u- (J) . 
e <* — 1 


t(cj+u/) 


t^cu+u/) 


(0 l 


(0l 


*u, in * tt.int 
a u ,L a u]',L 


~u,mt *u,m 
a u',L a u>,L 


|0i„) 


+ S (cj — U>') IOin) 


(53) 


Remarkably, the spectrum of quasiparticles emitted is precisely a bosonic thermal 
distribution with temperature a/ (27r); replacing fundamental constants, 


, m hot 

keT= 27- 


(54) 


The appearance of the 5 function shows that Eq. (53) represents a density rather 
than a number. This is because the orthonormal modes form a continuous spectrum, 
and are not normalized to unity but to a delta function - see Eqs. (41). An argument 


in pDj overcomes this problem by turning to wavepackets, and shows that the spectral 
flux density - the number of quasiparticles emitted per unit time per unit (angular) 
frequency interval - is obtained by dividing the calculated expectation value by 2i r: 

d 2 N 11 

(55) 


dtdu 


2vr e w - 1 


Due to the symmetry of the transformation (51), the expectation value of left- 
moving outgoing quasiparticles is exactly equal to the thermal spectrum of (53). 
There is a deeper significance to this than just symmetry, however. The trans¬ 
formation of the in- and out-operators allows us to write the in-vacuum explicitly 
in terms of the out-vacuum. Using the Fock basis, in which the annihilation and 
creation operators behave in the standard way [62] : 

Ou l">„ = Vn \n - 1) , aj, |n) = Vn + 1 \n + 1) , 


(56) 


the in-vacuum state is given by 
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where Z is a normalizing prefactor defined such that (0 
checked simply by acting on Eq 


|0 in ) = 1. (This can be 


(57) with an arbitrary in-annihilation operator, 
and utilizing the transformation (51), to show that it vanishes.) The fact that the 
out-creation operators appear only in R-L pairs shows that the radiation, though 
it looks thermal in the right and left sides separately, is strongly correlated between 
the two sides. Quasiparticles are emitted in pairs, one in the subsonic, the other 
in the supersonic region; and a measurement of the number of quasiparticles in any 
state on one side of the horizon infers that there are equally many quasiparticles in 
the corresponding state on the other side. The left and right systems, separated by 
the horizon, are maximally entangled [69] . This entanglement induces correlations 
between the left and right partners, which may prove invaluable in the eventual 
experimental detection of Hawking radiation [701 m] 72, E3j; see also <]T2.2| Most 


recently, the entropy associated with the entangled partners has been investigated 
in [731 EU - somewhat closing the circle of analogue systems by harking back to 
Bekenstein’s insight on gravitational black holes, but using the emitted radiation as 
a measure of entropy rather than the event horizon area. 


6 Discussion 


6.1 Summary 

Beginning with the generalized form of the black-hole spacetime (J7|) , we found the 
wave equation (12) for a massless scalar field, with the exact solution (flG| in terms 


of co- and counter-propagating components. The counter-propagating components 
were found to separate at the event horizon into two distinct spatial parts. Sta¬ 


tionary (single-frequency) modes were derived, satisying the dispersion relation (23) 
which combines the co-moving dispersion relation with the Doppler effect. An¬ 


other conserved quantity is the scalar product (25), which gives rise to positive- 


and negative-norm stationary modes related by complex conjugation. First assum¬ 
ing a homogeneous flow, the total field was written as a sum over the normalized 
stationary modes, and upon quantization it was found that positive- and negative- 
norm modes - corresponding to positive and negative co-moving frequency - are 
multiplied by bosonic annihilation and creation operators, respectively. Crucially, 
the relation between co-moving frequency and lab frequency changes sign between 
subsonic and supersonic flow, and while it is the former that determines the norm, 
it is the latter that is conserved. This leads, in the case of inhomogeneous flow with 
an event horizon, to the possibility of positive- and negative-norm mixing across the 
horizon. The bases of incoming and outgoing modes are generally different in any 
physical system, but the mixing of positive- and negative-norm modes leads natu¬ 


rally to a mixing of annihilation and creation operators (see Eq. (51)), and hence to 
the inequality of the in- and out-vacuum states which are defined by their respective 


annihilation operators (see Eq. (52)). So quasiparticles are spontaneously emitted, 


and are found to conform to the thermal spectrum (53), with temperature (54) pro¬ 


portional to the velocity gradient at the event horizon. Moreover, the quasiparticles 
emitted on each side of the horizon are maximally entangled, and emitted precisely 
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in pairs. 

The derivation has been kept quite general, in that it applies to any spacetime 
which can be described by the metric ([7]). Of course, this includes gravitational 
black holes: using V (r) = —Cyjrg/r, the velocity profile of the Schwarzschild black 
hole (see Eq. (§), we find that at the Schwarzschild radius V' (rg) = c/ (2 rg) = 
c 3 / (4 GM), and plugging into Eq. ( |54| ), we have k B T = he 3 / ( 8nGM) - exactly 
Hawking’s original formula, Eq. Q. 


6.2 Conceptual issues 


We identified the vacuum state as the in-vacuum, that which contains no ingoing 
quasiparticles; the ingoing quasiparticle states, in turn, were identified as the late¬ 
time equivalents of those corresponding to an initial horizonless spacetime. This 
inevitably raises an intriguing question: how is the vacuum state of an eternal 
horizon - with no horizonless spacetime in the asymptotic past - to be defined? 


There would seem to be no reason to use the modes (49) to define the vacuum state, 


since there can be no w-in modes in such a spacetime. However, the crucial property 
of the modes (49) is not so much their “ingoing” character, but their analyticity 
in the complex x-plane. One can be led to the same vacuum state by a simple 
and intuitive argument: that the state appears as vacuum to co-moving observers 
crossing the horizon [38] . This leads again to the requirement of analyticity across 
the horizon and hence to the “in”-modes (49) - a property selected, as remarked 


in 15.1, by the ingoing criterion in virtue of their analyticity in the horizonless 


asymptotic past. If the vacuum is defined differently, the modes used to define it 
would not be analytic at the horizon, and a freely-falling observer would observe an 
infinite flux of particles precisely at the horizon [38J. It thus appears that the state 
which gives rise to Hawking radiation is the only physically sensible one, and quite 
independent of the collapse phase. 

In either the collapsing or eternal spacetime, then, Hawking radiation seems 
to be intimately connected with the existence of wavepackets of the quantum field 
smoothly connecting the two regions which, in classical general relativity, are dis¬ 
connected by the event horizon. The horizon breaks these wavepackets into two 
pieces, which propagate out into their respective regions [^] Those pieces furthest 
from the horizon propagate outwards first, while those nearer the horizon escape 
later. As time progresses, then, the Hawking particles originate from wavepackets 
of ever-increasing wavevector, ever-increasing thinness, and take an ever-increasing 
time to propagate out from the horizon. In the dispersionless model, there exists 
an infinite reservoir of these trans-Planckian modes, so that the Hawking radiation 
can continue indefinitely (subject to eventual evaporation of the horizon, not ac¬ 
counted for here). It is a remarkably counter-intuitive feature of the dispersionless 
model that the late-time steady-state Hawking radiation originates from an infinity 
of modes “captured” at a single point of space (and time, if collapse is included), 
then gradually released and redshifted so as to produce a steady thermal spectrum. 


11 This is often interpreted as the creation of a particle-antiparticle pair via a quantum fluctuation, 
these then being separated by the horizon and unable to recombine mm- 
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The analytic continuation giving rise to the “in”-modes (49) is exact, even though 
it is derived from linearization of V at the horizon (performed in §3.3 to find the 
logarithmic phase). This linearization is easily justified: as discussed in §3.3 any 
wavepacket can be traced back in time to arbitrary thinness and arbitrary closeness 
to the horizon. The modes are thus well-defined arbitrarily close to the horizon, and 
the region over which analytic continuation onto the complex x-plane is imposed 
can be made arbitrarily small. Therefore, it does not matter how small is the region 
over which the linearization of the velocity profile is valid, so long as it is valid over 
some small region. Consequently, the Hawking temperature can depend only on the 
first derivative, a, of V precisely at the horizon. 

This picture breaks down entirely as soon as dispersion is introduced. If the 
wave speed is allowed to vary with wavelength, then the event horizon is not so 
well-defined as a single point and linearization of V at the horizon is not generally 
justified, so that we cannot hone in on a single parameter which uniquely determines 
the Hawking temperature. More drastically, a change of wave speed at some scale 
will tune wavepackets out of the grip of the horizon once they reach that scale. It 
is no longer true that any outgoing wavepacket must originate from one captured 
by the creation of the horizon; instead, all wavepackets, when traced back in time, 
will blueshift to the dispersive scale and then propagate away from the horizon, 
originating from a fixed finite wavevector at spatial infinity! Thus is the trans- 
Planckian problem resolved by dispersion - but can such different wave propagation 
yield similar effects to the trans-Planckian model? After all, given that a trans-sonic 
flow of dispersionless fluid can switch the sign of the norm by causing the wavevector 
to diverge at the horizon (see Figs. [I] and [5]), it is far from obvious that a similar kind 
of norm switching will still be possible in a dispersive fluid, and we might reasonably 
doubt the existence of Hawking radiation in dispersive systems at all. 
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Part II 

Dispersive model 


7 The dispersive wave equation 


Aiming to generalize the dispersion relation of waves in the medium, we are, due 
to the considerations of 


forced to abandon the spacetime metric as a starting 
point. We begin instead with the Lagrangian density, in which dispersion is modelled 
by the appearance of higher-order derivatives [211160]: 


£ = j l(9i + V<h)4>\ 2 - 1 |c (-id*) d x </>\ 2 


(58) 


The only difference from Eq. ( [ITj ) is the replacement of the constant c with the 
function c{k ); as a function of an operator, we consider it as a Taylor series in that 
operator. 


Application of the Principle of Least Action to Eq. (58) leads to the generalized 
Euler-Lagrange equation: 


dL d 


dL 


d 


dL 


d 2 


dt \d (d t (j)*) J dx \d(d x (j)*)J dx 2 \d (d%4>*) J dx 3 \d(d 3 (j)*) 


dL 


d 3 


dL 


Substituting L from Eq. (58), we find the wave equation EH 

( d t + d x V) (d t + Vd x )d~ c 2 ( -id x ) did - 0. 


+... = 0 . 

(59) 

(60) 


In the geometrical optics approximation, the velocity profile V (x) varies negligibly 
over a wavelength of the field q i, so that the “local wavevector” k(x) exists and 
is determined by the local value of V, analogously to Eqs. (22). This leads to a 


dispersion relation between u>, k and V: neglecting derivatives of V and k , we find 

(u - Vkf = c 2 (k) k 2 . (61) 


This is the generalization to dispersive media of Eq. (23). c(k) is the phase velocity 


LUcm/k in the co-moving frame, which describes the propagation of the phase kx — 
u>cmL in the lab frame, the phase velocity is rather uj/k. The group velocity is 
the derivative doj cm /dk = ±d(c(k)k)/dk in the co-moving frame, or du/dk = V ± 
d{c{k)k)/dk in the lab frame; this gives the velocity of the envelope of a wavepacket 
strongly peaked at the wavevector k. In the absence of dispersion, the phase and 
group velocities are identical. Typically, there exists a low-wavevector regime where 
the phase and group velocities can be approximated as constant, so that c ( k) —>• c 0 as 
k —y 0. For higher values of k, the dispersion can take two basic forms: superluminal 
and subluminal [j^] These are defined according to whether the magnitude of c (. k) 
becomes higher or lower than cq as k increases. 


12 Although these terms emphasise the Schwarzschild analogy, they have become standard terms 
to describe the dispersion in any context, much like the labels “subsonic” and “supersonic” to 
describe any type of flow. 


29 















As in §4.1 , there are two equivalent pictures for determining the possible solutions 
of the dispersion relation: from the point-of-view of the co-moving frame, they 
are the points of intersection of the dispersion profile ±\c(k)k\ with the co-moving 
frequency oj cm = oj — Vk', or, from the point-of-view of the lab frame, they are 
the points of intersection of the frequency u ; with the lab frame dispersion profile 
Vk± \c(k)k\. Also as before, there are two branches of the dispersion, and two ways 
to define them: 


• The u- AND u-BRANCHES: the sign of the phase velocity in the co-moving 
frame defines the counter-propagating w-branch (cu cm = c o — Vk = c(k)k) and 
the co-propagating t’-branch (ca cm = to — Vk = — c(k)k ); 

• The positive- and negative-norm branches: the sign of the frequency 
in the co-moving frame defines the positive-norm branch (cu cm = oj — Vk = 
\c(k)k\) and the negative-norm branch (a; cm = oj — Vk = —\c{k)k\). 


The definitions of the positive- and negative-norm branches carry over because, like 
its dispersionless counterpart the Lagrangian density ( [58]) remains invariant 
under phase rotations. Therefore, the scalar product of Eq. (25) is conserved even 
for dispersive media. 


8 Field decomposition with dispersion 


For a homogenous velocity profile, the wavevector A: is a conserved quantity, and the 
analysis of §4.3| — the definition of the orthonormal basis of fc-modes, the decomposi¬ 
tion of the held in this basis and the quantization of the held - carries over exactly, 
the only difference being the more complicated functional form of the frequency, 
uj(k). 

Frequency being the conserved quantity for a general inhomogeneous velocity 
prohle, we wish rather to express the held operator in the form of an integral over 


to: 


Hx,t) = 


dco 




(62) 


where the <p u (x) are the complete mode operators in the ^-representation, combin¬ 
ing all solutions of the same frequency. The transformation to the ^-representation 


performed in 14.4 was written for a general dispersion relation u{k), and therefore 
also carries over to the dispersive case. Dispersion complicates the transformation 
only through the introduction of additional counter-propagating solutions, which 
are easily found with the aid of dispersion diagrams. For simplicity, we shall re¬ 
strict ourselves to subluminal dispersion and examine its effects on the resulting 
transformation. Superluminal dispersion behaves analogously; see reference |59j for 
a detailed derivation of the transformation for both types of dispersion [^J Realistic 
dispersion relations may be complicated enough to introduce even more solutions 


13 Most of the analyses occurring in the literature (e.g., [701E33 Hi E3 EH) use superluminal 
dispersion - mainly because they are geared towards BECs, whose excitations follow a superluminal 
dispersion relation. 
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than are indicated hercj^J but careful attention to the dispersion diagrams should 
make the required modes clear. 

8.1 Homogeneous flow 

We learned from the dispersionless model that the set of plane-wave solutions to the 
dispersion relation depends crucially on the flow velocity - in particular, on whether 
it is subsonic or supersonic. Anticipating similar behaviour in the dispersive case, 
let us examine these two cases separately. 

Supersonic flow Suppose that the flow velocity is faster than Co, the low-wavevector 
limit of the phase velocity. The wavevector solutions of the dispersion are illustrated 
graphically in Figure [6] It is clear that, for all positive lab-frame frequencies co > 0, 
there exist only two wavevector solutions: a co-propagating t'-mode with positive 
norm and a counter-propagating u-mode with negative norm. This is exactly what 
was found for the same (supersonic) flow in the dispersionless model (see Figs. [4] and 
[5]), the only difference here being a different functional dependence of the wavevec- 
tors on frequency. The transformation to the ^-representation thus proceeds exactly 
as in §4.4] and the mode operators for oj > 0 are simply 

( x ) = ( x ) + (x) , (63) 

where the mode on the u -branch has been written as a Hermitian conjugate on ac¬ 
count of its negative norm. For oj < 0, the mode operators are simply the Hermitian 
conjugates of those in Eq. ( |63| : 0_ w (a;) = 4L{x). 

Subsonic flow Now consider a flow velocity slower than cq. The solutions are 
illustrated graphically in Figure [7J Here we observe behaviour significantly different 
from the dispersionless case. Firstly, let us note that, for all positive lab-frame 
frequencies 00 > 0, there exist a positive-norm u-rnode and a negative- norm w-mode, 
both being continuously connected to their counterparts in supersonic flow as V is 
varied from one value to the other. I 11 the dispersionless model, the negative-norm 
■u-mode experiences a diverging wavevector as the flow crosses from supersonic to 
subsonic (see Figs. [5] and [5]), becoming a positive-norm w-mode in the latter type 
of flow. Dispersion has completely removed this divergence, and consequently the 
negative-norm n-mode remains. 

However, we see from Fig. [7] that removal of the wavevector divergence does 
not quite prohibit the emergence of a positive- norm w-mode. More precisely, it is 
prohibited at high frequencies, above some threshold value cj max which depends on 
the value of the velocity (and vanishes in the limit where V becomes supersonic, 
so that the two cases are continuously connected). For frequencies below cu max , 
there comes into existence a pair of additional counter-propagating modes with 
positive norm. For convenience, let us label these ul and u2, where \k ul \ < \k u2 \; 

14 For instance, in the first paper [14] to consider dispersive analogue systems, liquid helium is 
considered as an example, with a dispersion profile that includes a roton dip. 
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the negative-norm w-mode, since it exists for all frequencies and velocities, will 
continue to be labelled simply u. The ul- and «2-waves can be considered as forming 
additional branches of solutions, of finite measure, which exist only at low frequencies 
and which vanish entirely when V becomes supersonic. Notice that, while the ul-, 
u2- and w-branches are all counter-propagating (i.e., they are right-moving in the 
co-moving frame), they are distinguished by their behaviour in the lab frame: 

• ul and u2 are distinct from u in that they have oppositely-signed lab-frame 
phase velocities oj/k ; and 

• ul is distinct from u2 in that they have oppositely-signed lab-frame group 
velocities doj/dk. 

The first property regarding the phase velocity also applies to the dispersionless 
case, since it distinguishes the ur- and UL-modes that exist on opposite sides of the 
horizon: it enables the mixing of waves of opposite norm. The second property, re¬ 
garding the group velocity, is a purely dispersive effect: it will be found to regularize 
the behaviour at the horizon. 

With these considerations, the mode operators can be written 


[x = 


a-X (®) + (a) + a u u4>f (x) + dl f Xl 0) for 0 < u < 


a u X (*) + a-Xl (x) 


for UJ > UJr, 


(64) 

Let us emphasise that the mode operators (64) for subsonic flow connect contin¬ 
uously to the mode operators (63) for supersonic flow, since cu max —>• 0 when V 
becomes supersonic. 


8.2 Inhomogeneous flow: Mode mixing 


As remarked in §4.4[ an inhomogeneous, asymptotically constant flow can be viewed 
as mixing or coupling the asymptotic plane waves to each other. Physically, this 
concept of mode mixing is perhaps better understood by considering the motion 
of a wavepacket, strongly peaked at a certain value of k. (Figures 11 and 12 in¬ 
clude results of numerical wavepacket propagation, alongside space-time diagrams, 
to help make the scattering concept more concrete.) In a region of constant V, the 
wavepacket will propagate in accordance with its group velocity, v g = dco/dk. A 
wavepacket incident on the inhomogeneous region will, after a complicated inter¬ 
action, scatter into some combination of resultant wavepackets, with frequency oj 
equal to that of the original wavepacket and wavevectors k corresponding to the 
possible solutions of the dispersion relation. To use the terminology introduced in 
§5.1[ there are two possible types of scattering mode: the out-mode, which contains 
a single outgoing wavevector in the asymptotic future (illustrated in Fig. 0 ; and 
the in-mode, which contains a single ingoing wavevector in the asymptotic past (il¬ 
lustrated in Fig. 12). Throughout the scattering process, the total norm is always 


conserved; but, as discussed in §5.2[ some of the waves involved may carry opposite 
signs of norm, leading to a degree of amplification. 
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Dispersive effects result in a richer variety of possible scattering processes, due 
to the appearance of additional plane wave solutions and hence more possibilities 
for mode coupling. Since the additional solutions exist only below the H-dependent 
critical frequency (u max , the scattering possibilities are themselves frequency- and 
velocity-dependent. The asymptotic flow velocities are generally different, yielding 
different critical frequencies for the left- and right-hand asymptotic regions. For 
definiteness, we denote these cu maX] i and where cu maX] i < (u maX!2 , and where 

one or both of the critical frequencies may be zero. These divide the spectrum into 
(at most) three regimes, in which quite different scattering behaviour is observed. 
Figures SI and [lO] show dispersion diagrams - with varying V to account for 
the inhomogeneous profile - for the various arrangements of the frequency regimes 
described below. Also note that V is assumed to be single-valued, so that only a 
single horizon (if any) exists. 


8.2.1 lo > U m M suppression of Hawking radiation 

This regime always exists for high enough frequencies, and in the case where the 
flow is everywhere supersonic so that both M ma r i and o ; maXi 2 vanish, it encompasses 
the entire spectrum; this situation is illustrated in Figure [8] 

In both asymptotic regions, the solutions of the dispersion are arranged similarly 
to Fig. [6] There exist only one positive-norm u-wave and one negative-norm u-wave. 
These vary continuously on their respective branches, being perfectly well-behaved 
at any horizon which might be present. This is quite different from the dispersionless 
case, where the u -wave experiencing divergence of its wavevector and a switch in its 
norm as the horizon is crossed. It appears there is no channel for Hawking radiation. 
It is not quite non-zero, though, since u-v mixing can occur, and since these have 
opposite norm they can form a Hawking pair [^J An interesting aspect of this pair 
is that, as can be seen in the dispersion diagram of Fig. [6| their group velocities 
are both negative, and both partners are emitted into the left-hand region. This 
is not the usual Hawking channel, though, and the strength of the u-v coupling 
is negligible for a slowly-varying velocity profile BED!- This high-frequency regime 
therefore corresponds to strong suppression of the Hawking radiation. 

Each of the u- and u-waves can form an ingoing or an outgoing wave. In either 
the in- or out-basis, then, there are two independent modes in this regime. 


8.2.2 o; maXi i < c o < LO the group-velocity horizon 

w m ax ,2 must be strictly non-zero for this regime to exist, so that at least one of the 
asymptotic velocities must be subsonic. If the other is supersonic, oo m * .v i vanishes 
and this regime extends down to zero frequency, a situation described by Fig. [9} 
On the other hand, if the flow is everywhere subsonic, cu maX) i is also non-zero, and 
this regime exists in a (possibly narrow) frequency window, as in Fig. 

In the asymptotic region corresponding to cu maXi2 , there are four wavevector so¬ 
lutions, while in the other asymptotic region there are only two. The two common 

lo This is not true for superluminal dispersion, where the u- and w-waves both have positive norm 

m- 
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solutions are the u- and u-waves; the extra solutions in one of the regions are nl 
and m2. Following the evolution of the dispersion diagram as V varies between the 
asymptotic velocities is particularly instructive; see Fig. |9j It is seen that the ad¬ 
ditional solutions nl and m2 vary towards each other, and, at some velocity Vg V h, 
merge into a single wavevector; as V is varied further, they cease to exist. Math¬ 
ematically, the ul- and u2-wavevectors cease to be real, but continue to exist as 
complex conjugates (as is most clearly seen in the complex fc-plane diagram of Fig. 
14); that is, they become evanescent waves - one exponentially damped, the other 
exponentially divergent - in the opposite asymptotic regime. Only the exponen¬ 
tially damped evanescent wave is allowed on physical grounds. Thus are we led to 
the conclusion that the ul- and u2-waves do not form two, but rather just one ad¬ 
ditional solution, for they must be combined in just the right combination to cancel 
the exponentially divergent wave [59]. We shall label this single additional solution 
nl2. 

The physical interpretation of this behaviour is again made clearer by consid¬ 
ering the evolution of wavepackets. In particular, we pay attention to their group 
velocities, v g = doo/dk = Vk± \c(k)k\, which is given pictorially by the slope of the 
lab-frame dispersion profile. As already noted, the ul- and n2-waves have oppositely- 
directed lab-frame group velocities: ul is right-moving in the lab frame, while m2 is 
left-moving. As the magnitude of V increases, the ul- and u2-wavevectors vary in 
such a way that the magnitudes of their group velocities decrease, until, precisely at 
the merging point, the group velocity vanishes. So a wavepacket on either the ul- or 
n2-branch, when sent towards the inhomogeneous region, will experience a shift in 
its wavevector such that it slows down to a standstill at the point where the velocity 
is Vg v h ■ It does not stay there, however, as numerical simulations show (see reference 
[ 21] and Fig. 0- Instead, its wavevector continues to evolve, crossing from one of 
these branches to the other, and its group velocity changes sign. (This is why they 
degenerate into a single branch: waves of one type evolve in time onto those of the 
other.) So the wavepacket continues to shift in wavevector, but moves back in the 
direction from which it came. The turning point where V = V gv } : is a group-velocity 
horizon, a dispersive analogue of the event horizon beyond which the wavepacket 
cannot propagate. It differs in that it does not cause the phase singularities of the 
dispersiouless model, and hence does not give rise to the trans-Planckian problem. 
But it differs also in that its position is frequency-dependent, so that it is impossible 
to speak unambiguously of “the horizon”. That said, in the limit oj —> 0 where the 
group velocity approaches Co, the group-velocity horizon also approaches a limiting 
point: the point where V = — cq . This low-frequency limiting horizon is the clos¬ 
est equivalent to the event horizon of the dispersionless model, since it ensures the 
existence of a horizon for all (counter-propagating) wavevectors in a spectral region 
around k = 0. 

Unfortunately, the simple geometrical picture just described does not capture 
the details of the full wave equation; in particular, it fails to capture the crucial 
coupling between the positive-norm nl2- and the negative-norm u-branches. Hawk¬ 
ing radiation can thus occur in n-nl2 pairs (as well as in u-v pairs). Whether the 
emitted quasiparticle has wavevector k ul or k u 2 depends on which of these is the 
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outgoing wave, which in turn depends on whether the velocity profile describes a 
black- or white-hole configuration: 

• Black HOLE: In this case, since V < 0, we must have Vr > Vl, and the 
ul- and «2-waves exist in the right-hand asymptotic region but not in the left. 
Then ul, with low wavevector, is the outgoing wave. The w-wave is always left- 
moving, and is thus emitted into the left-hand region. The Hawking partners 
are thus emitted in opposite directions from the group-velocity horizon, exactly 
as in the case of a gravitational black hole. 

• White hole: In this case, Vr < Vl, so the ul- and w2-waves exist in the 
left- hand region but not in the right. So it is u2, with high wavevector, which 
is outgoing. Moreover, it is emitted into the left-hand region along with the 
w-wave: the Hawking partners are emitted in the same direction! 

As earlier remarked, we restrict our attention mainly to black holes, but the white 
hole is included here for completeness, since its behaviour in this respect is so dif¬ 
ferent. Using super- rather than subluminal dispersion, however, does not yield 
drastically different results, merely switching the group velocities of the waves and 
hence the directions in which the quasiparticles are emitted [ 59] . 


8.2.3 0 < u) < ui. 


max,l • 


horizonless Hawking radiation 


If both asymptotic regions are subsonic (so that no low-frequency limiting hori¬ 
zon exists), then the critical frequencies are ordered thus: cu maXi2 > <Wiax,i > 0. 
The smaller critical frequency (u maXj i is strictly non-zero, and a third, low-frequency 
regime is available: 0 < oj < a wv ,i • It is shown in Fig. 10 
The ul- and w2-waves exist in both asymptotic regions 


| !85| . Following their 
variation with V, their wavevectors never reach the merging point; they experience 
no group-velocity horizon, and wavepackets in either of the two modes can propagate 
from one spatial infinity to the other. 

However, as we noted previously, the negative-norm u- wave exists regardless of 
the existence of the additional solutions. There is no need, as in the dispersionless 
case, to make a transition between subsonic and supersonic flow in order to have two 
solutions of opposite norm - they already exist, at a single value of V. So they can 
couple with each other to produce Hawking radiation, even in the complete absence 
of a group-velocity horizon [[85] ! 

Another novel aspect of this regime is that, there being no group-velocity horizon, 
ul and u2 do not degenerate into a single solution but remain as two completely 
separate solutions. They exist in both asymptotic regions, and so each of them 
can form an incoming wave or an outgoing wave. So each of them can couple 
to the negative-norm u-waves. There are thus three channels for the radiation, 
corresponding to the three distinct Hawking pairs u-v, u-ul and u-u2, each of which 
will be created according to its own spectrum. One of the latter two spectra will be 
a low-frequency continuation of that in the group-velocity horizon regime (ta maX! i < 
cu < cu maXj2 ), while the other exists only for oj < cu maXi i and must vanish outside 
this region. Generally, then, there are three Hawking spectra: the two (u max values 
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mark the cut-off frequencies of those involving only the counter-propagating waves, 
whereas the u-v pair has no cut-off. 


8.2.4 Form of the modes 

Taking the above considerations into account, the mode operators differ in the var¬ 
ious frequency regimes by the total number of independent solutions they contain: 



( (x) + 0) + a“ 2 0“ 2 (x) + (a) 

l (x) + d“ 12 0“ 12 (a) + (x ) 

[%,<% ( x) + (x) 


0 < u < (u maXi i 

w max ,i < u < u max,2 5 

CU > CU maX) 2 


(65) 

where we recall that u; maXi i or both u; max ,i and (n maXj 2 may vanish, depending on the 
nature of the asymptotic velocities. Note that the modes can be interpreted either 
as in-modes or out-modes, which exist in one-to-one correspondence. 


9 Hawking radiation in dispersive media 


9.1 Transforming between the in- and out-bases 


The (^-representation of the stationary modes given in Eq. (65) can be expressed as 
a dot product between a vector of operators and a vector of modes: 


iu (x) = (• • • 


/ : \ 

C’ in 

V : J 


( ■■ \ 


a: 


/x, out 


x,out 


V i 


( 66 ) 


where we have explicitly shown both the in- and out-mode representations. The 
generic modes and operators 0(j and stand for all the possible branches for any 
single value of uj, including any complex-conjugate modes and Hermitian-conjugate 
operators that may appear. (Recall that 0(j = 0^ and = d_ u ; these definitions 
take care of the negative-norm modes.) Since both the in- and out-modes form a 
complete set, they are related via a linear transformation: 
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( : 
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(67) 


Plugging this into Eq. 
tors: 


(66), we find a corresponding transformation for the opera- 


( : 

\ 


( : 

\ 



= S T 

~/x, out 
a u 


\ 

) 


V : 

/ 


( 68 ) 


where S T is the transpose of S. 
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S is simply the scattering matrix which describes how the outgoing modes scatter 
into ingoing modes when propagated backwards in time into the infinite past. Since 
the scalar product is bilinear and the in-modes are themselves orthogonal to each 
other, these coefficients can be written as a scalar product: S^ = (</£,, This also 
means that the inverse matrix - which describes the in-modes as linear combinations 
of out-modes - is easily obtained through use of the property (0 2 ,0i) = ( 0 i, 02 )*- 
In particular, the magnitudes of their elements are identical: S tw \ = [. In 

words: the absolute value of the amplitude of the u-in mode when the /i-out mode 
is scattered backwards in time is equal to the absolute value of the amplitude of the 
/i-out mode when the o-in mode is scattered forwards in time. This identity can be 


useful when calculating the radiation spectra, as we shall see in 111 


If some of the modes scatter into modes with oppositely-signed norm, then, 
when their amplitudes are substituted in Eq. (68), this will result in a mixing of 


annihilation and creation operators between the in- and out-bases. 


9.2 Spontaneous creation 


As in 15.2 we assume that the field is in the in-vacuum |0i n ) - that is, that there 


are no incoming particles, and the quantum state is annihilated by all in-mode 
annihilation operators. From Eq. (68), the expectation value of outgoing particles 
in a particular mode /a is 


<o,nl (asrfasy'io,.) = <o,„| 



|0,„) . (69) 


The subsequent algebra is entirely analogous to that in the derivation of Eq. (53). As 


there, the sum of operators that left-multiplies the in-vacuum state may be reduced 
to its creation operator terms, while the other is simply the Hermitian conjugate of 
this. Using the fact that the operators of different states commute, we have simply 


(0i„| ( 


«;°"') , < i £ ! ou ‘|0in) = 


Wl 




S (oji - u 2 ) 


(70) 


where u n denotes all states of oppositely-signed norm to the state /i. As before (see 


15.2), the occurrence of the <5 function shows that this is a number density. The 


spectral flux density - that is, the number of particles emitted per unit time per 
unit (angular) frequency interval - is obtained on division by 2tt 


d 2 N _ 1 
du dt 2n 




J'nM I 


(71) 


Recalling the definition of S in Eq. (67), the spectral flux density of particles in the 
mode /i is proportional to 

• the sum of the squared amplitudes of the opposite-norm in-modes into which 
the /i-out mode scatters when propagated backwards in time; or 
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the sum of the squared amplitudes of the //-out mode over all opposite-norm 
in-modes when they are scattered forwards in time. 


Equation (71) is the generalization to dispersive media of Eq. (53). There are 


three important differences to note. Firstly, due to the inability to factorize the 


dispersive wave equation (60) into u- and v-parts as in Eq. (15), it is now possible 


for u- and v-waves to couple to each other. This coupling is negligible in the limit of a 
slowly-varying velocity profile [40], but in general it should be taken into account (see 
references (mug and [50J for predicted u-v coupling in Bose-Einstein condensates). 
Secondly, pairs created in dispersive media are not necessarily localised to opposite 
sides of the horizon. In the usual black hole case, it is certainly true, for the created 
ul- and u -waves have opposite group velocities. But for a white hole [53], or in the 
horizonless case where the ul- and //2-waves exist in either asymptotic region [85], it 
is possible to create a u-u2 pair, with the same sign of the group velocity and hence 
emitted into the same asymptotic region^} Finally, unlike the dispersionless case, 
the scattering amplitudes in Eq. (71) are generally unknown, and must be found 
either approximately or numerically. 


10 Calculation methods 

The dispersive framework just developed does not lead to an exact analytical deriva¬ 
tion of the Hawking spectrum, as in Eq. (53). There are no exact solutions akin 


to (16), valid for arbitrary velocity profiles. The main exact solutions we do have 


are the plane waves (44) whenever V is constant; Ending how these are coupled by 


inhomogeneities in V requires either numerical or approximate methods. 


There are two arbitrary functions in the wave equation (60): the velocity profile 
V ( x ) and the dispersion profile c(k). These are defined in the dual spaces of x and fc; 
therefore, we should be mindful that, depending on their functional forms, it might 
be more convenient to work in one space than the other. If one of these functions 
appears in the wave equation as a finite polynomial, for example, then it will form 
an exact finite-order differential equation in the dual space. Since V(x) is taken to 
be asymptotically constant, only c(k) can be described as such, leading to an exact 
differential equation in x-space which can be solved numerically. 

On the other hand, if V (a;) is monotonic, then the appearance of multiple so¬ 
lutions means that the position of a given wavevector is better-defined than the 
wavevector at a given position. Moreover, the group-velocity horizon is also better- 
defined in terms of positions of wavevectors, which vary quite smoothly in contrast 
to the abrupt transition from real to complex wavevectors that the horizon engen¬ 
ders in rc-space. From an analytical point of view, then, the problem is better-suited 
to the k -space representation. 


16 Studying the emission from a white hole in BEC, in [53] it was found that interference between 
the correlated u- and u2-quasiparticles leads to a very pronounced “checkerboard” pattern in the 
density-density correlations. 
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10.1 Numerical method: differential equation in x-space 


We shall first consider the intuitive method of solving the wave equation in position- 
space through numerical integration, thus yielding exact solutions. This method is 
also used in EflEaiZaEDlEZI. 

Assuming a stationary solution (p UJ (x)e ~ lut , the wave equation (60) is transformed 
into its time-independent form 


(-•w + d x V) (-iw + Vd x ) <t> w - e 2 (—id x ) dl<j> = 0 . 


(72) 


This is a linear ordinary differential equation. If the phase velocity c 2 (k) forms a 
finite polynomial in k, then the differential equation is of finite degree and may 
be integrated numerically provided appropriate boundary conditions are specified. 
Different boundary conditions lead to the various solutions associated with a given 
frequency. 

As always, V is taken to be asymptotically constant. In the asymptotic regions, 


then, the normalized plane waves (44) are particular solutions, and the general 


solution is a linear combination of these: 


( x ) = c^ 


x 


(73) 


where p is a generic label for the various plane-wave solutions, and the c M are arbi¬ 
trary complex constants. Some of the may have negative norm, and can always 
be written as the complex conjugate to emphasize this. Others may have com¬ 
plex wavevectors, and since the dispersion relation is real, these always occur in 
complex-conjugate pairs. Only that wavevector decreasing towards infinity is phys¬ 
ically allowed, so only it can appear with non-zero coefficient in the general solution 


(73). Note that, since complex wavevectors do not correspond to propagating waves, 
their “norm” is physically meaningless (it is actually zero |5lj). and their amplitudes 
do not enter into equations of norm conservation. 


The general solution (73) applies in each of the asymptotic regions separately , 


with modes and coefficients particular to their own region. The left- and right-hand 
solutions are not independent, but linearly related through the linear differential 
equation ( |72| . We may think of the set of plane-wave solutions in one region as a 
basis of solutions, and the velocity profile as a linear operator, transforming between 
the left- and right-bases. This is entirely analogous to the linear transformation 
between the in- and out-bases, Eq. (67), except that the spatial form of the left- 


and right-modes makes them better-suited to solution by numerical integration. 
Specifying boundary conditions appropriate to a single plane wave in one asymptotic 
region, Eq. (72) can be numerically integrated through to the other region, and its 


solution there fitted to the sum of plane waves (73); this is illustrated pictorially in 
Fig. 13 Working through all possible initial plane waves, the transformation matrix 
- a transfer matrix , linking two regions in space - is built up column by column. 
This yields the matrix T such that 
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while for the coefficients we have 


/ = \ 


0 /x,left 
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= 7 ^ 


r,v, right 


i / 


(75) 


where 7^ is the transpose of T. 

Any in- or out-mode can be solved for in terms of its scattering coefficients using 
Eq. (75). We simply set to zero all those coefficients corresponding to ingoing or 
outgoing plane waves other than the one we are considering, as well as the coefficients 
of any exponentially diverging solutions [^J There remains the coefficient of the 
single incoming or outgoing wave, which is set to unity to normalize the complete 
mode. Plugging these constraints into Eq. ( |75| yields the scattering coefficients of 
the mode. The norm of each individual scattered wavevector is simply the squared 
modulus of its scattering coefficient, with the inclusion of a minus sign for those 
with negative norm. Conservation of norm implies that the sum of the norms of the 
scattered waves should equal unity, the norm of the single ingoing or outgoing wave; 
this is very useful as a numerical check. 


10.2 Analytic method: the step-discontinuous limit 


While numerical integration as just described can be applied to an arbitrary velocity 
profile, there exist two limiting regimes in which instructive approximation methods 
can be applied. These regimes relate to the steepness of the velocity profile in relation 
to the fundamental length scale determined by the dispersion. If the velocity changes 
over a distance much longer than this dispersive length scale, then the system is in 
the slowly-varying regime of geometrical optics; an analytical method applicable to 
this regime is described in detail in §10.3 On the other hand, if the velocity changes 
over a distance significantly shorter than the dispersive length scale, the waves of 
the system are unable to resolve the details of the change, so that shortening the 
transition region even further will have no significant effect. In this regime, the 
system is well-approximated by a velocity profile with a step discontinuity [^] This 
type of profile is studied in references [.751 IS ES (SBj. 

The step-discontinuous velocity profile is exactly solved by a sum ( |73j ) of plane 
waves (and non-divergent evanescent waves) in each of its constant-velocity regions, 


1 ' If the divergence of a complex wavevector is too strong, it may swamp the numerics and make 
the calculated solution unreliable [59j. The solution can be tested by checking whether the norm 
is conserved, and whether the solution is significantly affected by increasing the accuracy and 
precision of the numerical integration. 

18 A detailed study in [£jj[] shows that the relevant dispersive scale to be used in such a comparison 
is actually the cut-off frequency w max , and not - as might be expected - the wavevector at which 
the dispersion profile deviates from the dispersionless model. 
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and these solutions must be matched appropriately at each discontinuity. The 
matching conditions are determined by the differential equation (72): as the dis¬ 


continuous limit is approached by an ever-steepening profile, Eq. (72) must be 


satisfied throughout the steepening process. So long as the system is dispersive so 


that Eq. (72) contains spatial derivatives of order greater than 2, this limit is well- 
defined. Suppose for definiteness that the differential equation is of order 2 n, the 
phase velocity being described by the order-(2n — 2) polynomial 


n—1 


c 2 (k)=J2 c jk 21 


(76) 


3 = 0 


Then the stationary field 0 and its first 2 n — 2 spatial derivatives are continuous 
everywhere, while at a point of discontinuity where V changes abruptly from V~ to 
V + , the (2 n — l)-th derivative of 0 has the discontinuity 

(9j”“V) + - (9« 2 ""Vr = T^~ M - < V+ + V-)d^} (V + - V~) . (77) 

'-'n— 1 

Since 0 and d x (f) are assumed continuous everywhere, there is no ambiguity in their 
values at the point of discontinuity. These conditions - the continuity of 0 through 
<9 2n ~ 2 0 and the discontinuity condition (77) for <9 2n_1 0 - provide 2 n independent 


relations between the 2 n unknown amplitudes of the ingoing or outgoing waves (rel¬ 
ative to the amplitude of the single outgoing or ingoing wave, typically set to unity). 
This is then a solvable linear system. Since there is no need to perform numerical 
integration or to build up the transfer matrix T by considering one asymptotic plane 
wave at a time, the solution can be found much faster than in the general case of 


110.1 


Condition (77) is particular to the system obeying wave equation (60), though 
other systems will admit analogous discontinuity conditions E What is noteworthy 
is that, in the presence of dispersion, this infinite-steepness limit exists and can be 
calculated, for the reason outlined above: dispersion introduces a fundamental length 
scale beneath which the details of the velocity transition cannot be resolved. The 
Hawking flux must therefore approach a finite limit in the limit of infinite steepness. 
This is in stark contrast to the dispersionless case, where the Hawking temperature 
is simply proportional to the steepness at the horizon, and would thus diverge in 
the discontinuous limit. From a physical viewpoint, the absence of dispersion allows 
infinite resolution, so that a decrease of the length of the transition region will always 
be “visible” to the waves of the system, and will therefore affect the Hawking flux. 
This is yet another rewording of the trans-Planckian problem, and the existence 
of the infinite-steepness limit is another manifestation of the regularization of this 
problem brought about by dispersion. 


19 Reference ESI gives a similar such condition in a one-dimensional Bose-Einstein condensate, 
while works out condition ( |77| ) in detail and gives a similar derivation for nonlinear light 
interaction in optical fibres. 
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10.3 Analytical method: the slowly-varying limit 

The earliest analytical descriptions of Hawking radiation in dispersive media [22} fT8 j 
tackle the problem in fc-space, where horizons are better-defined, and transform 
back to :r-space via the saddle-point approximation. The solution in fc-space is 
approximated by linearizing the velocity profile in the vicinity of the horizon. This 
is a severe limitation of the method, for even if the geometrical optics approximation 
holds, it is not necessarily true that higher-order derivatives at the horizon are 
negligible. Other treatments also rely on linearization at the horizon [79, [80], or on 
a specific form for V [TO]. 

Here, an alternative theory is presented (as also in references [BZ} Elj). Inspired 
by earlier treatments, it is developed mainly in fc-space, and uses the saddle-point 
approximation to transform back to x-space. But V is not linearized; it is completely 
general, except that it is assumed to be slowly-varying in such a way that a WKB- 
type approximation is justified. 


10.3.1 Solution in fc-space 


We take the Fourier transform of Eq. 
x —* id k and d x —* ik. Rearranging sligh 
/c-space takes the form 


(|72|) by making the substitutions 0 —>• 0, 
ly, we find that the wave equation in 


OJ 


j, - V(id k ) 


c\k) 


keb ) = 0 . 


(78) 


The velocity profile can be characterized by a scale, and we indicate this by including 
a parameter e in its argument: V(id k ) —» V{ied k )\ this is later included in the 
definition of V by formally setting e — 1. The purpose of e is as a measure of the 
steepness of V, so that a “slowly-varying” velocity profile can be defined as that 
remaining in the limit of small ePl We make the ansatz k<b = exp (0), where bp may 
be Laurent-expanded in powers of e: bp — e -1 (y 0 + e 7h + + •••)• Plugging this 


ansatz into the left-hand side of Eq. (78) yields a Taylor series in e: 


OJ 


-~V{ied k )) — c 2 {k) 


+ e 


OJ 


--V{id k ip 0 )) — c 2 {k) 


d k <Pi + ^d k In (V'{id k tp 0 )c{k)) 


+ 


kcj) 


Equation (78) is solved by setting the coefficient of each power of e to zero, yielding 
an infinite sequence of coupled equations for the coefficients bpj. We implement the 
slowly-varying approximation for V by truncating the Taylor series at the first power 


20 In a mathematical sense, e plays here a role analogous to that played by h in the WKB 
approximation in Quantum Mechanics. Just as the WKB approximation works well in the semi- 
classical limit, the present approximation works well in the geometrical optics limit; both cases 
refer to a “slowly-varying” background. 
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of e. This leaves only the two lowest-order equations: 


(^-V(id k tpo)) — c 2 (k) = 0, 

(79) 

d k Vi + ^d k In ( V\id k yo)c{k )) = 0 . 

(80) 


Equation (79) is precisely the dispersion relation (61), except that the flow ve¬ 
locity V is not assumed to be constant. The argument of V is the position x, but 
imposing Eq. (79) selects a particular position for each value of k: that position 
at which A: is a solution of the dispersion relation for the local value of V. Let us 
denote this wavevector-dependent position by x(k), and note its distinctness from 
the position variable x. Explicitly, x(k) is defined such that 


(r - V(x(k)j) 2 - <?(k) = 0, (81) 

or, inverting this relation, 

x (k) = V^ 1 - c(k)^ , (82) 


Note that we have chosen the sign of the square root of c 2 (k ) to correspond with the 
counter-propagating n-modes; changing this sign gives another, independent, func¬ 
tion corresponding to the w-modes, but we shall not consider this here 21 Equation 


(79) is solved by setting idx-fo equal to x(k), or, integrating, 


<A) = -i / x(k')dk' + C 0 , 


(83) 


where the lower limit of integration is omitted as it can be incorporated into the 
constant of integration Cq. 


The second of our sequence of equations, Eq. (80), may now be solved by sub¬ 
stituting x(k) for idkfio, and it is easily seen that 


^--Miawow^ + c,. 


(84) 


Plugging the solutions for tp 0 and <pi into the original ansatz, we find the approximate 
k -space solution 


0 ( k ) « 


kJV'( X (k))c(kC P 




X(k')dk' , 


(85) 


where the constants of integration Cq and C\ are incorporated into the overall pref¬ 
actor <f> 0 . 

21 The independence of the position functions corresponding to the u- and ti-modes shows that, 
to this level of approximation, these modes are decoupled from each other. So, while we could 
perform the following analysis for n-modes, we would not predict any scattering from u- into v- 
modes or vice versa. This corroborates the fact that the u-v coupling becomes negligible in the 
slowly-varying regime [40j . 
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The economical value of the fc-space representation is demonstrated by Eq. (85). 


Whereas, in the x-space representation, the number of independent solutions can 
(depending on the dispersion relation) be arbitrarily large, in A;-space we have found 
only two independent solutions, corresponding to the counter- and co-propagating 
modes; neglecting the latter, the solution is uniquely determined up to an unim¬ 
portant multiplicative prefactor. The question inevitably arises: how can these 
dramatically different dimensionalities be reconciled, when the two representations 
ostensibly describe one and the same system? The answer is that, in Fourier trans¬ 


forming back to x-space, the fc-space solution (85) can yield different x-space so¬ 


lutions depending on the integration contour taken in the complex /e-plane 
Just as the precise linear combination of x-space solutions is selected by appropriate 
boundary conditions, so the integration contour in /c-space is also selected by appro¬ 
priate boundary conditions; but the latter - especially in the presence of a horizon 
- is much easier to identify. 


10.3.2 Transforming to x-space: Saddle-point approximation 


Inverse Fourier transforming Eq. (85) to find the x-space solution, 


(p(x) = 


r+o o 

' — oo 
/• + 00 


(t> (k) e lkx dk 
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Loo k \/V'(x{k))c(k) 

f -|-oo 


exp ^ixk — i J x(k') dk^j dk 


A(k) exp (iO(k)) dk , 


( 86 ) 


where we have defined 
A(k) = 


•hr 


WV'(x(k))c(k) 


, 6[k) = xk 


X(k') dk'. 


(87) 


Its integrand taking the form of an amplitude multiplied by a phase factor, the 
integral in Eq. (86) seems to lend itself well to the saddle-point approximation. 


This observation is corroborated by noticing that, keeping the scale parameter e, it 
would appear only as a factor of 1/e multiplying the phase 9{k) (with x rescaled 
by e so that ex remains constant). So, in the limit of small e - the limit in which 
our approximation for <fi(k) is valid - the saddle-point approximation also becomes 
a valid means of evaluating the inverse Fourier transform in the asymptotic regions. 

The saddle points of the phase 6{k ) are simply those points in the complex k- 


plane at which the derivative dkO(k ) vanishes. From the second of Eqs. (87), we see 
that such points occur where x — y(fc) = 0; recalling the definition of \(/c) in Eqs. 


(81) and (82), this condition simply inverts the relationship between position and 


wavevector so that, for a given x, the saddle points lie at those values of k which are 
solutions of the dispersion relation for the local value of V. As x is varied from one 
asymptotic region to the other, the saddle points trace out all values of k which solve 


the dispersion relation (81) for real x\ in particular, they vary continuously between 
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wavevector solutions in the left-hand region and those in the right, at which extremes 
x(k) is singular. We shall denote the various wavevector solutions by fiT( x ), so as 
not to confuse them with the wavevector variable k. 

Consider the example of Figure 14, which shows the imaginary part of \(k) 
for the example velocity and dispersion profiles specified later in { 11.1 The lo¬ 


cus of saddle points (equivalently, of wavevector solutions) for real x is shown as 
blue dashed lines. This forms two disconnected structures. On the right, with 
real positive k, are the positive-norm solutions, intersected by complex-conjugate 
wavevector solutions. This crossbow-like structure is symptomatic of the presence 
of a group-velocity horizon: at the point of intersection, where the group velocity 
vanishes, two real wavevector solutions merge and form complex-conjugate wavevec¬ 
tor solutions beyond the horizon. The real wavevectors on either side of the group- 
velocity-horizon point are precisely the ul- and w2-branches. On the left, with real 
negative k, lies the w-branch, which does not experience a horizon and therefore 
exists throughout the real space, in agreement with §8.1[ 

The saddle-point approximation treats the integral of an amplitude times a phase 
factor - as in Eq. (86) — as a sum of contributions from the saddle points of the 
phase. In the vicinity of a saddle point, the phase is Taylor-expanded to second 
order: 

1 , , , ,. / ,,2 


6(k) 


8(^(x)) + ^dld(^ l (x)) (.k - ^(x))' 


( 88 ) 


Substituting back into the integrand of Eq. (86), and assuming that the amplitude 
varies slowly enough so that it can be replaced by the constant value A(k^(x)), the 
integrand takes the simple Gaussian form: 


A(k) exp (■ iO{k )) « A{k^{x)) exp (id(K^(x))) exp (i^dl0(K^(x)) (.k — x )) 2 

(89) 

Upon integration, this is easily evaluated; but before we can do so, we must check 
that the resulting integral is convergent. It will only be so if the integration contour 
can be deformed to pass through the saddle point in its direction of steepest descent, 
i.e., the direction in which the exponent i\d^.6{n^{x)) (k — k^(x )) 2 is negative. In 
this direction, the integrand is indeed a convergent Gaussian function, at least in the 
vicinity of the saddle point; and, so long as the magnitude of the phase increases fast 
enough (equivalently, so long as e is small enough, or the velocity profile is sufficiently 
slowly-varying), then the Gaussian function converges to zero quickly so that only 
the behaviour in the vicinity of the saddle point is important in the evaluation of 
the integral. (Retrospectively, this justifies replacing the amplitude A(k) with its 
value A(k,^(x)) at the saddle point.) Then the integral over this relatively small 
region can be treated as an integral over the entire real line, the differences far from 
the saddle point having no significant effect. So a single saddle point, for which the 
contour can be deformed to lie along the direction of steepest descent, contributes 
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approximately 


*4.(/ffi(x)) exp (^(^(x))) 


r >+00 


exp ( i-dlO(K fl (x)) ( k — x )) 2 j dk 


= A(k^(x)) exp (i 9 (k^{x))) 


2vr 


i dl9(n f) '(x)) 


( 90 ) 


to the integral of Eq. (86). If the integration contour passes through several such 


saddle points, these contributions are added together. 


From the definition of \{k) in Eqs. (81) and (82), and after some algebraic 
manipulation, the second derivative of the phase is given by 


d 2 k d(^(x)) = -d k x(^(x)) = 


Vg(X, X )) 


(91) 


K»(x)V'(x) 

Here, v g (x,k ) = V(x)+ d k {c{k)k) is the group velocity of the wavevector k at the 


position x. Applying directly to Eq. (86), substituting the saddle point contributions 


(90) and replacing the amplitudes and phases of Eqs. (87), we find 


<t>{x) 




2vr 


i dl9(n^(x)) 


exp (i 9(n^(x))) 


<h 0 V2ttY1 


±i- 

e 4 


-> ^ ^ exp / 


^^(^(x)) K,^{x) Vg(x, ^(x))] 

/t^(±(X>) 


exp i x k^(x) — i 




7r 


x(k ) dk ± i— 



In the last line, the solution has been restricted to the asymptotic constant-velocity 
regions, x —> ±oo. There, the solutions fiT become independent of x, and are equal 
to the wavevector solutions for the given asymptotic flow velocities. The various 
plane wave solutions and their fc-dependent amplitudes are found to be, up to a 
constant prefactor, the u;-normalized plane waves of Eq. (44). So the solution in 


the asymptotic regions is, as one would expect, just a linear combination of these 
various plane wave solutions. What is remarkable about Eq. (92), though, is that 


the relative coefficients of these plane wave solutions are not arbitrary: they are 
determined by the phase integral of x{k), which in turn encodes the effect of the full 
velocity profile on the scattering between the different wavevector solutions. The 
real part of this integral gives the relative phases between the plane waves; and, since 
x(k) can be analytically continued into the complex fc-plane, the integral can also 
contain an imaginary part, yielding their relative amplitudes. The only ambiguity is 
in the path taken in the complex plane between the various wavevectors, or, more 
specifically, its topology with respect to the singularities and branch cuts of x(fc). 
This is exactly as noted previously: in the fc-space representation, the selection of 
a potentially large number of coefficients of x-space solutions is replaced by the 
selection of a single contour in the complex /c-plane. 
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10.3.3 Selection of the contour 


As previously stated, the integration contour is selected so as to satisfy appropriate 
boundary conditions. These boundary conditions relate to the “in” or “out” nature 
of the mode we are aiming for, for this determines that certain plane waves - those in- 
or outgoing waves which do not correspond to the single in- or outgoing wave of the 
mode - must have zero amplitude. We must therefore single out the particular mode 
that relates to Hawking radiation. Neglecting the u-branch, the outgoing waves are 
ul and u, so the radiation is produced in ul-u pairs and the relevant modes are 
the w 1-out and w-out modes. Of these, the former is the easiest to calculate, for in 
the presence of a group-velocity horizon and for both subluminal and superluminal 
dispersion, it has a purely evanescent character in one of the asymptotic regions. 
(This mode is illustrated by a space-time diagram and a wavepacket simulation in 
Fig. 11 ) The wl-out mode thus contains only an exponentially damped solution 
beyond the horizon. This is the boundary condition which selects the required 
integration contour. 

Let us return to Fig. 14 Firstly, we note that, to ensure convergence, the 
integrand of Eq. (86) must tend to zero at the limits of the integral. This restricts 
the possible directions in which the integration contour can tend to infinity. The 
simplest (though by no means the only) way in which this condition can be satisfied 
is to have the contour tend from — oo, where the imaginary part of x(/c) is positive, 
to Too, where the imaginary part of x(fc) is negative. In Fig. [l4j we see then that 
there is a tendency of the contour to lie on the lower-half of the complex fc-plane. 

Now consider the boundary condition of pure exponential damping. The wavevec- 
tors of interest are those in the asymptotic regions, which lie at the extremities of the 
locus of saddle points shown as blue dashed lines. These have been colour-coded, 
with red representing the wavevector solutions in the right-hand subsonic region, 
and black representing the wavevector solutions in the left-hand supersonic region. 
It is in the left-hand region, beyond the horizon, that the evanescent modes exist. 
There, the exponentially damped solution has negative imaginary part: it is that 
complex wavevector lying on the lower-half complex plane, at the lower end of the 
crossbow. Also shown in Fig. 14 are the directions of steepest descent through the 
extreme saddle points 22 For the exponentially damped wavevector, this direction 
is almost horizontal. The contour, then, is that from — oo to Too which lies entirely 


on the lower-half plane. It is shown as a white curve in Fig. [14] 

The specification of the contour determines the complete solution in x-space 
through insertion into Eq. (92), because the contour can be continuously deformed 
so long as its topology with respect to the singularities and branch cuts of x(&) is 
not changed. That is, we may move the contour in the complex fc-plane, but so 
long as we do not cross any of the branch cuts where x(^) is discontinuous (where 
light shading meets dark shading in Fig. 14), the result of Eqs. (86) and ((921) do 


not change. The final check, then, on whether the selected contour is the correct 
one is to see whether it can be deformed so as to pass through the other (colour 


22 The directions of steepest descent are actually calculated in the limit as x —> Too, for the 
extreme points themselves are singular and the directions of steepest descent are ill-defined precisely 
at those points. 
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black) saddle points in their directions of steepest descent. It cannot; therefore, only 
the exponentially damped wavevector contributes to the solution in the left-hand 
supersonic region, and the selected contour is indeed that which corresponds to the 
wl-out mode. 


10.3.4 Extracting the Hawking temperature 


Recalling Eq. (70), the Hawking flux is determined by the squared amplitude of 
the oppositely-normed ingoing wave relative to that of the outgoing wave. For the 
wl-out mode, as shown in Fig. 11, the scattering process takes place entirely in 
the right-hand subsonic region. So, having selected the integration contour using 
boundary conditions at x —* — oo, we now consider the solution at x —> +oo. The 
saddle points now lie at the wavevector solutions u, ul and u2 coloured red in Fig. 


14 


Which of these wavevectors contribute to the solution is determined by whether 
the integration contour can be deformed to pass through these saddle points in 
their directions of steepest descent; it is easily seen that it can be deformed to pass 
through all three of them. The x-space solution in this region is therefore a sum of 
the three plane waves, one outgoing (ul) and two ingoing {u2 and u), with relative 


amplitudes given by Eq. (92). The significance of the integration contour is not yet 


spent, for it determines the path taken between the three wavevector solutions in 
the complex fc-plane when calculating the phase integral of x(k). 

Denoting the amplitudes of the ingoing waves u2 and u as a u and /3 W , respectively, 
and assuming that the overall mode is normalized so that the outgoing wl-wave has 
unit norm, we have, by norm conservation, 


cu 


-iA.r=i. 


(93) 


The Hawking flux is determined by the amplitude j3 u of the u-wave, the negative- 
norm ingoing component. Defining the frequency-dependent temperature T(oS) in 


accordance with the form of the Planck spectrum of Eq. (53), 


w = 


ghw/ksT^u)) 


(94) 


and utilizing Eq. (93), we have 


ol. 


Pu 


= exp 


hu 


k B T(uj) 


(95) 


,/& 


the amplitude of the «2-wave relative to the w-wave - is, according 

23 1 of the phase integral of x(k) from k\ to 

In 


Now, | 

to Eq. ( 92h , given by the imaginary part 


/ch 2 along a path consistent with the integration contour selected in f 10.3.3 


2 i 


23 The real part of the integral gives their relative phase, which is unimportant in calculating the 
flux. 

24 Note that the integral taken between k p 1 and k^ is purely real, so that, according to this 
prediction, the incoming u2-mode and the outgoing ul-mode have equal amplitude, and the ex¬ 
ponentiated imaginary part of the integral taken between the positive- and negative-fc branches 
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Fig. 14, such a path lies in the lower-half complex /e-plane, an important detail in 


that it determines how the central branch cut is to be circumvented. 

Although the result could be left formulated in this way, it is made more elegant 
by noticing that the imaginary part of the phase integral can be isolated by closing 
the contour in the upper -half complex plane. Being an analytic continuation of 
an essentially real-valued quantity, the wavevector-dependent position y(fc) must 
obey the relation \{k*) = So, if we were to choose the contour to lie on 

the upper-half plane - equivalently, transforming the integration variable k to its 
complex conjugate k* - this would simply result in the complex conjugate of the 
original integral. But, on closing the contour in this fashion, we subtract this value 
from the original integral, so that the real part is cancelled out while the imaginary 
part is doubled. Since the phase integral appears as an exponent in Eq. ( [92] ) , this 
closed phase integral gives precisely the squared relative amplitude la^/^l 2 , and by 


comparison with Eq. (95) we have 


Tiu) 


k B T(u) 


X(k)dk 


(96) 


where the closed integration contour is taken around the central branch cut, as 
shown by the black curve in Fig. 


14 


Equation (96) is a generalization of the dispersionless result (54), which, as we 
now show (see also [8T]). it reproduces. The phase integral of Eq. (96), instead 


of an integral of position over wavevector, may also be written as an integral of 
wavevector over position by integration by parts: § x(k)dk — — <f> n(x)dx. While 
in the general dispersive case the branch points and the integration contour take 
complicated forms in position space (see, e.g., references |5T1 f54] for a position- 
space analysis of Bogoliubov modes in BEC), the dispersionless case is very simple 
there because it has only one singular point: the event horizon itself, where the 
wavevector diverges. This is easily seen by examining the closed form of the position- 


dependent wavevector (recall Eq. (22)): k u (x) — oj/(V(x) + c), which is analytic in 


the complex x-plane except at the point where V(x) = —c, where it has a simple 
pole. The closed contour integral is equal to 2tt times the residue of this pole; writing 
V{x) = —c + ax + ..., it is seen that the residue is co/a. Then we have 


hu 

k^r 


n u {x)dx 


—dx 
ax 


2ttuj 


a 


(97) 


from which we retrieve Eq. (54). 


Equation (97) also holds in the low-frequency limit of the dispersive model, for 


then the central branch cut closes in closer to the origin, and since the dispersion 
profile is assumed to approach c(k) = Cq in this limit, we can then treat the phase 
integral as if the system were dispersionless. So, as u —> 0, Eq. (96) reproduces 


can equally well be interpreted as \a u //3 u \ or as simply |/3 W |. This is a defect of the approach, 
and we observe that of these two possibilities, the phase integral should give \a ul /(3 ul \, since this 


relates the two high-magnitude wavevectors where the amplitude A(k) (defined in Eqs. (87)) is 
better-behaved. Other approaches combining WKB techniques in momentum- and position-space 
are able to resolve this issue - see references m , m and [S3]. 
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precisely the dispersionless result (54), where the horizon is understood as the low- 
frequency limiting horizon where V(x) = — Cq. This can be taken as a measure of 
the applicability of Eq. (96): if the low-frequency temperature T(u — 0) is equal to 
that predicted by the dispersionless model, then the phase integral of Eq. (96) can 
be used to find the temperature for higher frequencies. We shall examine a model 
in q fTT| in which this is found to be the case, and we shall see how the applicability 
of Eq. (96) is related to the criterion of a “slowly-varying” velocity profile. 


11 Application to a simple model 

Apart from restricting ourselves to dispersion relations of the subluminal type, the 
dispersive framework of previous sections has been kept fairly general, and we have 
noted only the main differences occurring when superluminal dispersion is used. 
In this section, let us turn to a specific example, providing a quantitative account 
of, but mainly to get a feeling for the qualitative nature of the expected Hawking 
radiation. The velocity and dispersion profiles used are kept as simple as possible, 
so that the qualitative properties discovered may still be expected to apply quite 
generally. 


11.1 Specifying the model 

11.1.1 Dispersion profile 

The dispersionless case corresponds to a wavelength-independent phase velocity: 
c 2 {k ) = Cq. Since dispersive effects arise through the occurrence of higher-order 
derivatives in the wave equation - and bearing in mind that c 2 (—k ) = c 2 (k ) for an 
isotropic medium - the simplest deviation from the dispersionless model is through 
the inclusion of a quadratic term in c 2 (k) [^J 

c 2 (k) = c 2 0 ^l - . (98) 


This dispersion relation is illustrated in Figure |f5[q). Note that it reduces to the 
dispersionless form when k 2 jk\ <C 1. The parameter kd, then, indicates how large 
k must be before dispersive effects begin to be felt. 

In line with the framework we have so far considered, we have chosen the disper¬ 


sion relation (98) to be sub luminal: the phase velocity of short wavelengths drops 


below the long-wavelength limit Co- It is also possible to include the dispersive term 
with a plus sign, in which case the dispersion would be superluminal. Indeed, a 


25 Dispersion relation (98) - also used in [60] - has a strict cut-off wavevector kd , at which 
the free-fall frequency vanishes. This differs from the subluminal dispersion relations shown in 
Figs. m where a limiting free-fall frequency is approached at high k , similarly to the dispersion 
considered in m■ So long as the flow velocity does not become too small - as will be the case 
here - there is no significant difference in the behaviour of these two forms, since the wavevector 
remains well within the cut-off. For low asymptotic velocity, the wavevector may reach kd-, at 
which point this dispersion relation breaks down [60] • However, even the flat dispersion relation 
becomes problematic in the low-velocity limit; see “Conceptual issues” of (12.2 and §VI of [BO]. 
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quadratic superluminal dispersion relation - exactly of the type (98) with the sign 
switched - is obeyed by the excitations of Bose-Einstein condensates [51] . and has 
been numerically studied quite extensively in j59j EH 85; [86] . We shall find that the 


results obtained from subluminal dispersion are qualitatively very similar. 


11.1.2 Velocity profile 

Endeavouring to keep the physical system as simple as possible, we shall use a 
velocity profile that varies monotonically between two asymptotic values; such a 


profile shall then be invertible, and can be used in the WKB-type analysis of £ 10.3 
We choose 

V (x) = ^ (V n + V L ) + ^ (V R - V L ) tanh {ax) . (99) 

This is illustrated in Figure [l5](6) . Vr and Vl are the asymptotic values of V ( x ) as x 
tends to +oo (the right-hand region) or —oo (the left-hand region), respectively. The 
parameter a characterizes the steepness of the velocity profile, but this also depends 
on the difference Vr — Vl', a more direct interpretation of a is that it governs the 
length of the region of transition between Vl and Vr, which is of order 2/a as shown 
in Figure [l5|(b) . As a —> oo, the velocity approaches the discontinuous step-profile 


v^oo( x) = v L e{-x) + v R 0{x ), 


( 100 ) 


where 9{x) is the Heaviside step function, equal to zero for x < 0 and unity for 
x > 0. In this limiting case, the treatment described in §10.2 may be applied. 

We restrict our attention to a left-moving and accelerating flow, i.e. Vl < Vr < 0. 
Therefore, if a horizon is present, it corresponds to a black hole horizon. White hole 
spectra are studied in 


11.1.3 Normalizing the wave equation 


With the dispersion and velocity profiles given by Eqs. (98) and (99), the wave 


equation (60) contains five parameters: c and kd from the dispersion profile; and Vl, 
Vr and a from the velocity profile. However, two of these correspond to the scaling 
of space and time; in redefining the variables to make them dimensionless, we can 
reduce the number of independent parameters to three. 

Let us proceed by defining a dimensionless velocity, which is the flow velocity 
divided by the low-wavevector phase velocity: U = V/cq, and also by defining a 
dimensionless distance, which is the “phase” of a wave with wavevector kd : X = 
kdX. Then a natural unit of time is l/c^kd, and we define the dimensionless time 
T = Cokdt. Dimensionless wavevectors and frequencies are defined accordingly: 
K = k/kd and = uj/c 0 kd . With these definitions, phases, which are anyway 
dimensionless quantities, are unchanged: kx = KX and out = QT. In terms of the 


dimensionless variables, the wave equation (60) becomes 


{dr + dxU) {8 t + U dx) 4> — C 2 (—idx) dx4> — 0 j 
where the rescaled dispersion relation is 

C 2 ( K ) = 1 - K‘ 2 . 


( 101 ) 


( 102 ) 
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and the rescaled velocity profile is 

U (X) = \{Ur + U l ) + \{U r - U l ) tanh (aX) . (103) 

The wave equation now contains only three parameters, and these relate only to the 
velocity profile U(X): 

• Ur = Vr/cq and Ul = Vl/cq, the asymptotic values of the flow velocity relative 
to the low-wavevector wave speed; and 


a = a/kd, a parameter that combines steepness of the flow velocity profile 
with the wavevector characterising dispersion. Since the length of the region 
over which V changes is on the order of 2/ct, we see that 1/a is approximately 
the phase difference of the wavevector kd across this transition region. So, if a 
is small, this transition region is many cycles of kd wide; whereas, if a is large, 
V varies quickly over only a few cycles of kd- The geometrical optics and step 
function limits of § §10.3 and |10.2] correspond to the limits of small and large a, 
respectively; the transition between these regimes will be discussed in §11.2.3 
(see also 


Although in principle the values of these parameters can be set arbitrarily, there 
will typically be practical limits. This is particularly relevant for the steepness 
parameter a, since the same dispersive effects which regulate the singular nature of 
the horizon will also tend to smooth out the sharpness of the velocity transition. 
In BEC, for example, they are regulated by the so-called healing length, which 
determines both the scale of short distance dispersion and the minimum length of 
the transition region, a — 1, then, represents a value which would be difficult to 
exceed in practical 

The difference in the asymptotic velocities can, in some contexts, also be limited. 
Nonlinear optics provides a good example: a light pulse induces an effective velocity 
difference on the order of the change in refractive index, typically no more than 
10~ 3 [SSI [57]. There would seem to be no generally applicable limits to the velocity 
values, however. 


11.1.4 Required modes 


Before we can use the 5-matrix of (|9]to calculate the quasiparticle creation rates, we 
must first select the appropriate modes whose scattering amplitudes encode their 
values. For the case of a black-hole velocity profile with group-velocity horizon, 
we have actually already found such a mode in §10.3.3 since it is the boundary 
conditions appropriate to the desired mode which there determine the integration 


contour in momentum space. This was the rt 1-out mode - illustrated in Fig. 11 
- which happened to have a single exponentially damped solution in the left-hand 
supersonic region. Now, whereas the analytical approximation of [ 10.3 treats u- and 


rvmodes as decoupled, this is not true in general, and through numerical integration 


26 In principle, however, arbitrarily high steepness can be achieved by suitably combining an 
external potential with a spatially-varying interaction constant mmm- 
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the coupling is included automatically. We should thus also account for the u-rnode. 
As can be seen in Fig. [lTl the ingoing u-wave appears in the right-hand subsonic 


region along with all the other propagating waves. So it remains true that the full ul- 
out mode is exactly equal to the left-mode which is purely a damped exponential in 
the left-hand region 2 ‘, and it may be found by straightforward integration, without 


any subsequent linear transformation. Since the u -wave is the only one with opposite 
norm to ul, the creation rate of the ul-u pair is (as found in 19.2) given by the 


relative squared amplitude of the ingoing w-wave compared to that of the outgoing 
wl-wave. 

There is another possible creation channel for this setup: the u- and u-waves 
have opposite norm, so u-v pairs can be created. Their creation rate is given by 
the relative squared amplitude of the ingoing u -wave to that of the outgoing u-wave 
in the u-out mode, ffowever, since the u-wave is always left-moving, the outgoing 
t?-wave is in the left-hand supersonic region, where the evanescent waves also exist. 
There can be a problem with numerical stability when we try to isolate the left- 
hand u-wave, since numerical integration in the left-hand region - no matter what 
the initial conditions - will tend to lead to exponential growth [Bill .59] • That said, 
in the limiting step function profile there is no such problem in calculating the 
u-out mode, since this involves only matching at a single point and no numerical 
integration. 

Finally, consider the case when there is no group-velocity horizon. Then all three 
of the waves ul, u2 and v can be outgoing, so all three of them can be created by 
pairing with the u- wave. As we have just been considering, their creation rates 
could be calculated by calculating each of their out-modes and singling out the 
relative squared amplitude of the ingoing u- wave, ffowever, recall from §[9] that 
this is equivalent to summing their relative squared amplitudes over all in-modes 
of negative norm. This latter approach is the most economical because there is 
only one such in-mode, corresponding to the u-wave. It is shown in Fig. 12, where 


it can be seen that it contains only the u- and wl-waves in the right-hand region, 
and no ul-wave in the left-hand region. Two integrations are required to calculate 
it: we find first the til-right and tj-right modes, determine the amplitudes of their 
left-hand components, then form the linear combination which exactly cancels the 
til-wave in the left-hand region. From the resulting mode, the creation rates are 
the squared amplitudes of the outgoing ul-, u2- and u-waves relative to that of the 
ingoing u-wave. 


11.1.5 Predictions 

Before turning to numerical results, let us first extract the analytical predictions of 


the techniques of 110.2 and 110.3 


In the high-steepness limit where the velocity profile approaches that of Eq. 


(100), we simply need to solve a linear system involving derivatives of the plane 
waves, and hence powers of K for each possible wavevector solution. Since the 


2 'This is not true for superluminal dispersion, where the ingoing u-wave appears in the same 
region as the evanescent wave 155] , 
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dispersion (Q — UK) 2 = C 2 (K)K 2 is of fourth order, there will be four such solutions 
and the linear system is 4 x 4. The resulting expressions are unwieldy; but, in the 
low-frequency limit where the Planck spectrum becomes the simple pole T/Q p*j 
the expression for the u-ul Hawking spectrum also approaches this form, and the 
low-frequency temperature can be read off. After some straightforward but tedious 
algebra, we find 


ri-ihigh a _ tj2\ 

1 fl^0 ~ V 1 U R) 


1/2 (1 + U R ) (-1 - U L ) (~Ur ~ Ul) 
(1 - Ur) (1 - U L ) (U R - U L ) 


(104) 


A more detailed derivation can be 


This result was first stated in reference 
found in Appendix A of reference [87]. 

In the low-steepness limit where V can be treated as slowly-varying, the tem¬ 
perature at any frequency should be well-approximated by the phase integral of Eq. 


(96). Moreover, this phase integral takes a very simple analytic form when the ve¬ 
locity profile is of the hyperbolic tangent form (99), as was used to plot Fig. [I4| at 


the branch cuts, the imaginary part of x(k) jumps between — 7t/2 a and +7r/2a, while 
the central branch cut around which the integral is taken lies exactly on the real 
axis. We can deform the contour to lie exactly along the branch cut, in which case 
it becomes simply the real integral of a constant. After some trivial rearrangement, 
Eq. (96) then gives (in units where Ti — k B — 1) 


rplow a /ry\ _ 

K) -*K#-Kt 


it 


(105) 


where K^ and Kf are the (normalized) ul- and w-wavevectors calculated in the 
right- and left-hand regions, respectively. It is easily shown (see §4.6 of reference 
[87] ) that, in the low-frequency limit, this becomes 


a (1 + U R ) (-1 - Ul) 
" tt (U r - U l ) 


(106) 


which, as can be checked by evaluating the derivative of V at the horizon, is exactly 
Hawking’s prediction, Eq. (54). As observed at the end of 110.3, Eq. (96) always 
yields this result in the low-frequency limit. 


11.2 Calculated spectra 


The richness of the dispersive spectra is perhaps best appreciated by fixing one of 
the asymptotic velocities and varying the other through U = — 1. This allows us 
to examine spectra from the various frequency regimes described in 
manner in which they vary between these regimes. 


.2 and the 
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The normalized frequency-dependent temperature is defined by analogy with Eq. (53) such 


that the relative squared amplitude of the corresponding mode is given by the expression for the 
Planck spectrum: |dn| 2 = (exp(f 1/T(Q,)) — l) -1 . 
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11.2.1 Varying Ul'- the subsonic side 


In Figure 16, the parameter a is fixed at 0.1, Ur is fixed at 
from —0.51 to 


-1.5; Figure 17 


-0.5 and Ul is varied 
on the other hand, shows spectra for the same 
values of Ul and Ur but in the discontinuous limit a —> oo. (We call this the 
“subsonic side” because the central value of the velocity is greater than —1.) When 
Ul is close to Ur, the spectral region that experiences a group-velocity horizon (i.e., 
cU mav i < co < cu max , 2 ) is narrow, and most of the spectrum is in the horizonless 
(c o < (U max i regime. As Ul drifts away from Ur, cu maXj i decreases so the spectral 
region experiencing a group-velocity horizon widens at the expense of the horizonless 
regime; the group-velocity horizon regime encompasses the entire spectrum when Ul 
reaches —1 and cu maXj i vanishes. As Ul varies further, we are in a true “black hole” 
regime, where a low-frequency limiting horizon (V = — Co) exists. 

In Figures 16 and 17 the solid lines plot the frequency-dependent temperature 
of the wl-wave, while the dashed lines do so for the w2-wave, which is only emitted 
where there is no group-velocity horizon (i.e., where 17 < 17 max ,i)- The dotted lines 
of Figure [T7| plot the temperature of the r?-wave, which is difficult to calculate for 
finite a = 0.1 and, where it is possible, is found to be vanishingly small. 

A list of the noteworthy features of the spectra follows. For the ul-waves: 


• The ul spectra show a general increase in temperature as Ul is decreased 
(i.e., as the left-hand flow is made faster). Since the derivative of the velocity 
profile increases in this direction, this is to be expected. As the limit of 
infinite steepness, the discontinuous limit shows temperatures higher by about 
an order of magnitude compared with the finite value a = 0.1. 


• The low-frequency temperature of the ul spectra increases in the same fashion, 
but it is finite only when the low-frequency limiting horizon exists (i.e., when 
Ul < —1)- On the other hand, when the flow is everywhere subsonic, the 
low-frequency temperature vanishes. This is the main effect of the presence 
or otherwise of the horizon; in particular, the spectrum at higher frequencies 
varies quite smoothly over Ul = — 1. 


• The spectra typically show significant dependence of temperature on frequency, 
and are therefore far from Planckian. Interestingly, the most Planckian spec¬ 
trum (i.e., that with the flattest temperature curve) occurs for Ul = —1.5, 
where the velocity profile is symmetric around the horizon Ul = —1- Indeed, 
the more asymmetric the velocity profile, the greater the deviations from the 
Planck spectrum seem to be. 

• In the low-steepness regime and when Ul > — 1 (so that U is everywhere 
subsonic), the ul spectra vary smoothly between the horizonless and group- 
velocity horizon regimes (i.e., at o; max ,i)- In the discontinuous limit, however, 
the transition between these regimes at o; ma . v i is more pronounced. 

• There is a cut-off in the ul spectra at cu maXi2 , at which point the outgoing 
ul-waves cease to exist. In the low-steepness regime, this cut-off is very sharp; 
in the discontinuous limit, it becomes more rounded, the fall-off beginning at 
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lower frequencies. Since the wl-waves are emitted into the right-hand region, 
their cut-off frequency depends only on Ur, and since this is held fixed all ul 
spectra are seen to vanish at the same frequency. 

For the w2-waves: 

• At a fixed frequency, the u2 spectra also show increasing temperature with 
decreasing Ul, and temperatures higher by about an order of magnitude in 
the discontinuous limit. Since the u2 spectra exist only when there is no 
low-frequency limiting horizon, their temperature always vanishes in the low- 
frequency limit. 

• The u2 spectra experience a cut-off at cu maXj i, at which point the outgoing 
w2-wave ceases to exist and a group-velocity horizon comes into existence. 

For the u-waves: 

• The v spectra, shown only in the discontinuous case, generally have lower tem¬ 
peratures than the ul- and w2-waves (for equal frequency and Ul)- Following 
the general trend, their temperature increases with decreasing Ul, though only 
significantly so when Uj j is close to Ur. Their low-frequency temperature al¬ 
ways vanishes, in accordance with the fact that the u-waves never experience 
a horizon. 

• Since the u- and u-waves exist even when the ul- and w2-waves do not, the 
r^-spectra do not experience any cut-off. 


11.2.2 Varying Ur: the supersonic side 

Consider now the case where Ul is fixed at a supersonic value and Ur is varied 
through —1. (The central velocity will always be less than —1, hence we call this the 
“supersonic side”.) Since the ul- and w2-waves do not exist when the flow velocity 
is supersonic, there can be no mixing and hence no radiation when Ur < —1, and 
the spectrum is expected to vanish in this limit. Moreover, as Ur decreases (i.e., the 
flow in the right-hand region becomes faster), the cut-off frequency of the wl-waves 
decreases. Therefore, the spectrum should become narrower with decreasing Ur, 
approaching zero width in the limit Ur —> —1 where the wl-waves cease to exist. 

In Figure [18J a is fixed at 0.1, Ul is fixed at —1.5 and Ur is varied from —0.5 to 
—0.9; in Figure [l9j Ul and Ur take the same values while the spectra are calculated 
in the discontinuous limit a —> oo. There is always a point at which U = —1, so 
there is no outgoing w2-wave and only the ul and v spectra are non-zero. 

These spectra behave exactly as predicted: they both become narrower and 
experience an overall decrease in temperature as Ur is decreased. Notice that, as 
before, the discontinuous limit shows temperatures about an order of magnitude 
greater than for a = 0.1, and that the spectra become very rounded so that the 
fall-off takes hold at lower frequencies than it does for a = 0.1. 


In the discontinuous limit (Fig. 19), we also include the spectra of the u-waves. 


As noted previously, these have generally lower temperatures than the wl-waves, 
with zero temperature in the low-frequency limit, and do not experience a cut-off. 
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11.2.3 Varying a: the low- and high-steepness regimes 


Finally, let us examine the transition between the low- and high-steepness regimes 
by continuously varying a. Rather than consider the evolution of a whole spectrum, 
let us consider the evolution of the low-frequency temperature in the limit O —y 0. 
As evidenced in Figs. [16] and 17, for this to be finite requires the existence of a 
low-frequency horizon where U = — 1, and so Ul < — 1 and Ur > — 1. 


Figure 20 plots the low-frequency temperature against a/ao, where ao is the 


value of a at which the low-steepness prediction, Eq. (|106|), is exactly equal to the 
high-steepness prediction, Eq. 


(104): 


a 0 = 7T (l - Ur) 


1/2 {-U R - U L ) 

(1 - U R ) (1 - U L ) 


(107) 


The high- and low-steepness predictions are also shown in Fig. (20), and by con 


struction they meet at a/ao = 1- It can been seen that a/ao provides a very good 
indicator for the regimes of low steepness (where the low-frequency temperature 


obeys Hawking’s formula and the spectrum as a whole obeys Eq. (96)) and of high 


steepness (where the step-function approximation of (10.2 is valid) 


11.2.4 Comparison with superluminal dispersion 

It is instructive to compare these observations with what is reported in the literature 
for superluminal dispersior 29 The most relevant paper is [85], since it also considers 


highly asymmetric profiles and includes the regime of all-supersonic flow with no low- 
frequency limiting horizon (although it considers fixed velocity difference Ur — Ul 
and a change of the central velocity, so the actual temperature values are not directly 
comparable). The resulting ul spectra are qualitatively very similar to those of Fig. 


16[ when a low-frequency limiting horizon where U — — 1 is present, the spectra 
approach finite temperatures in the low-frequency limit O —^ 0; and all spectra show 
a cut-off at some maximum frequency. An important difference due to the reversed 
behaviour of the w-modes for sub- and superluminal dispersion with respect to sub- 
and supersonic flow (see ([8]) is in the reversed behaviour with respect to asymmetry 
towards the sub- and supersonic sides: 

• When the velocity profile is asymmetric and towards the supersonic side, the 
spectrum shows strong enhancement towards the high-frequency end. This 


is analogous to what we have found for subluminal dispersion in Fig. 16 for 
velocity profiles asymmetric towards the subsonic side. 

When the velocity profile is asymmetric and towards the subsonic side, the 
spectrum narrows and is suppressed at higher frequencies, just as we have 


seen in Fig. 18 for velocity profiles asymmetric towards the supersonic side. 


29 The reader is reminded that subluminal and superluminal dispersion for the wave equation (601 
are studied and compared in (551. Also, §IV E of [831 details an approximate symmetry between 
subluminal and superluminal dispersion, helping to explain many of the observations made here. 
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Reference [85] does not calculate temperatures for the case of no group-velocity 
horizon, but it does find non-vanishing spectra for the case of an entirely supersonic 
flow in the mid-frequency regime where a group-velocity horizon exists, just as we 
did in Fig. 16 for the case of an entirely subsonic flow. 

The v spectra are not calculated in |85]. However, we can predict that these will 
behave very differently in the cases of subluminal and superluminal dispersion: 


• For subluminal dispersion, the v-wave has opposite norm to the u- wave; so, it 
can form a Hawking pair with the u -wave in the same way that the u 1- and 
w2-waves can, and since there is no cut-off frequency for the mnodes, this v 
spectrum exists in principle for all frequencies. 

• For superluminal dispersion, the r;-wave has opposite norm to the ul- and u2- 
waves; it thus forms Hawking pairs with these in the same way as the w-wave, 
and since the ul- and u2-waves have cut-off frequencies above which they no 
longer exist, so these v spectra will vanish at these cut-off frequencies. 


The transition between the low- and high-steepness regimes was recently studied 
in detail in [86] in the context of Bose-Einstein condensates. A smooth transition 
between the two regimes is also observed there. It is found that, to determine in 
which regime the system lies, one should compare the steepness of the transition 
(the “surface gravity” n) with o; max , the cut-off frequency of the spectrum. This 
suggests that co max /k is simply related to a/a 0 , the normalized steepness used to 
plot Fig. 20. Work on this particular aspect is still ongoing, however. 


12 Discussion 


12.1 Summary 


Generalizing the Lagrangian (10) derived from the spacetime metric to the La- 


grangian (58) exhibiting dispersion, we find a wave equation (60) and a dispersion 


relation (61) analogous to the dispersionless model but with a fc-dependent phase 


velocity c(k ), which we take to approach cq as k —» 0 and which might be superlumi¬ 


nal or subluminal at higher wavevectors. The scalar product (25) is also applicable 


with dispersion, and on writing the total held as an integral over stationary modes 
and quantizing, it is still true that positive- and negative-norm modes are multiplied 
by bosonic annihilation and creation operators, respectively. However, the switch in 
the norm of the counter-propagating mode in the dispersionless model is replaced 
by a different kind of behaviour. There is no divergence of k anywhere, so there 
is a counter-propagating wave which exists for all values of V, never changing the 
sign of its norm. However, for either subsonic or supersonic how (depending on the 
type of dispersion), there comes into existence a pair of counter-propagating waves 
of opposite norm to the already existing one. This allows opposite-norm pairing 
as before, though - thanks to shifting and reflection at a group-velocity horizon 


without trans-Planckian frequencies in the past (though see fl2.2 for a caveat). 


More surprisingly, it also allows opposite-norm pairing in the complete absence of a 
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group-velocity horizon, and since this can occur in two possible pairings, there are 
two types of Hawking spectra in such a situation. The spectra can be calculated 
numerically by solving for the appropriate out- or in-modes via numerical integra¬ 
tion; or, when the steepness is sufficiently large, by solving a linear system with 
matching conditions at a discontinuity; or, when there is a group-velocity horizon 
and the velocity gradient is not too great, by a phase integral in fc-space. Assuming 
a simple model of fluid flow and dispersion, we have found such Hawking spectra, 
confirming their theoretical existence and that they agree with the phase-integral 
formula in the appropriate regime. 

12.2 Overview of current research 

The field of analogue Hawking radiation has recently experienced its own inflationary 
period, the number of interested groups and proposed experimental setups having 
increased dramatically. Here 1 offer a brief overview of the current efforts and 
challenges in this field. For a review and complete historical overview of the general 
field of analogue gravity, see [ IS] . 

Analogue systems The earliest system to be considered for analogue Hawking 
radiation was superfluid helium nn 122], being the main experimentally-accessiblc 
liquid in which quantum effects are important [^j Once dilute gas Bose-Einstein 
condensates became experimentally accessible, their use as analogue systems was 
soon realised [23112l] : their theoretical simplicity pT] and the relative ease with which 
they can be controlled in laboratory conditions have made them the most widely 
studied of all the analogue systems. That said, due to flow instabilities induced by 
supersonic flow, it is not so easy to actually create a black hole configuration in 
superfluid systems, and it was only recently that this was finally achieved in BEC 

m- 

A change of reference frame can give rise to a black hole analogue in a station¬ 
ary medium, so long as there is an effective “velocity profile” which is moving with 
respect to the medium. This was first considered in the context of solitons in super¬ 
fluid helium ESI, and has also been considered for solitons in BECs [92]. It forms 
the basis of the analogue model in nonlinear optics: a strong pulse of light induces 
an effective change in the refractive index of the material, and this perturbed index 
travels through the medium with the pulse [551 931191) 95] • In the co-moving frame 
of the pulse, the effective index is stationary while the material is moving - exactly 
the type of situation studied in this tutorial, with the labels “co-moving” and “lab” 
switched. This can give rise to a group-velocity horizon, and the resultant shift 
between the ul- and w2-waves has been measured in optical fibres [88]. Advantages 
of these systems are the high effective velocity - which Doppler shifts the emitted 
photons into the UV range, making them easier to detect - and the self-steepening 
effect, which can form optical shocks with very high steepness [551 57]. Their main 

30 Superfluid helium-3 provides an entire array of cosmological analogies nanni, in which ana¬ 
logues of the field theories familiar to us appear as low-frequency limits of more fundamental 
theories. 
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disadvantage is the smallness of the effective velocity change induced by the pulse 
[87 , 153] , Given the complicated dispersion relations of optical media, the details 
of the Hawking process in these systems exhibit a rich phenomenology, as recently 
studied in [95] . 

Although they are far from the quantum regime with no prospect of emitting 
detectable Hawking radiation, analogue black holes have been formed by trans-sonic 
flows of water gZl SSI S3!, and the scattering behaviour of incident waves has been 
observed. In particular, the stimulated Hawking effect - the classical scattering 
between positive- and negative-norm waves - can and has been observed in this 
system [47, 39]. 

We have taken only a glimpse at the wealth of proposed analogue systems, which 
also includes: ions moving on a ring mm, “photon fluids” in microcavities [98] 
and in an analogue Laval nozzle [99] : and “fluids” of exciton-polaritons [ 1001 101J . 


Detection of Hawking radiation Creating an analogue system is one thing; 
detecting Hawking radiation is quite another, and has so far proved elusive. The 
predicted Hawking temperature of the setup realised in reference [91] is on the orders 
of nanokelvin - on the same order, but still somewhat less, than the temperature of 
the condensate itself. The Hawking spectrum will then be practically indistinguish¬ 
able from the thermal background. This problem is reinforced when losses are taken 
into account [102 , 103], with quasiparticles induced by losses typically overwhelming 
the Hawking quasiparticles [1031]. 

A promising proposal to overcome this thermal swamping is to exploit the entan¬ 
glement between the Hawking partners, which leads to correlations between them. 
These correlations have been studied extensively in BECs [S3 EH E21E3], where 
they show up as clear patterns in density-density correlations. Remarkably, the cor¬ 
relation signatures remain strong even after the imposition of a finite condensate 
temperature higher than the Hawking temperature. 

Another technique to increase the visibility of the Hawking spectrum is to make 
it highly non-thermal by inducing sharp resonance peaks. This requires a several- 
horizon setup so that interference between them occurs. If the asymptotic regions are 
one subsonic and one supersonic - as studied in the BEG context in [ 104 ] - the system 
is (usually) dynamically stable, so that direct observation of spontaneous Hawking 
radiation is possible, its non-thermal character making it easier to distinguish from 
actual thermal radiation. The black hole laser (previously mentioned in [2.4) - 
of experimental relevance since many practical setups include a black-hole white- 
hole pair - also yields strong resonant peaks at well-defined frequencies [551 EB* 52, 
158] : however, this radiation is strongly coherent and essentially classical in nature, 
forming from a runaway amplification process of stimulated emission (as in a laser) 
that drowns out any spontaneous component of the radiation. 

A signal claimed to represent Hawking radiation has been measured emitted from 
laser filaments in bulk nonlinear media, at a wavelength conforming to theoretical 
predictions [ 105] . The interpretation of the signal, however, is controversial, with 
claims that the results are not consistent with Hawking radiation and that they 
are signs of another kind of creation process related to the non-stationarity of the 
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background |106l 1107] . 


Theoretical methods Similarly to the theoretical results presented in §11] there 
are a number of works which investigate Hawking radiation by mode analysis in order 
to infer the spectrum for given dispersion and velocity profiles; references mm 
provide the earliest examples. Recent avenues of investigation include: comparison 
of black and white hole radiation for sub- and superluminal dispersion [59]; spectrum 
and density-density correlations from a step-discontinuity in BEC [76]; the influence 
of asymmetric velocity profiles [85]; the dependence of the spectrum on multiple 
length scales m, the transition between the low- and high-steepness regimes [86]: 
and mode analyses of black hole lasers [57] [58] - 

In a manner similar to that of §10.3[ there are several works which attempt purely 
analytic derivations of the Hawking spectrum in dispersive system. The earliest 
examples are [221 GE], which recover the Hawking result to leading order. This 
was subsequently generalized: to general preferred reference frames ESI; to high- 
frequency corrections [80]; to more general dispersion relations [82] : to correlation 


functions [83]; and, as presented in [10.3, to more general velocity profiles [81j- The 
latter approach was also considered in [87] for the “horizonless” case, where y(/c) 
has two central branch cuts in the complex fc-plane rather than just one as in Fig. 
14 there are signs that this case does obey some kind of generalized phase-integral 
method, possibly involving linear combinations of various contours. 

Analytical methods are almost always approximate to some degree, while mode 
analyses focus on stationary states (unless unstable modes are included, as for black 
hole lasers). In order both to check our analytical results and to gain insight into 
the time-dependence of the system, ab initio simulations are invaluable. These 
have been used in the context of BECs, in the evaluation of the density-density 
correlation pattern [ZDEIJ and on the dynamics of the periodic configuration with 
a black hole-white hole pair [54]; the growth of an oscillating “checkerboard” pattern 
in the white hole case [53] is an interesting example of an unexpected feature revealed 
by simulations. 


Conceptual issues It was early recognised [3HHB0] that, in the case of a Schwarzschild 
black hole (or indeed any velocity profile with vanishing asymptotic velocity), dis¬ 
persion does not quite solve the trans-Planckian problem. Though it regularizes 
the wavevector divergence at the horizon by the shifting of high-wavevector ingoing 
waves to low-wavevector outgoing waves, tracing the ingoing waves back in time to 
spatial infinity sends them to arbitrarily high wavevectors (MIEDJ- So the position 
of the divergence has been shifted to infinity, but the divergence still exists in prin¬ 
ciple. A very closely related point is that of conservation of the lab frequency when 
V —>■ 0 at spatial infinity, since then the lab and co-moving frames coincide in the 
asymptotic region. Thus, conservation of the lab frequency is equivalent there to 
conservation of co-moving frequency; and, since the sign of the co-moving frequency 
gives the sign of the norm (see Eq. (29)), it follows there can be no conversion 
between positive- and negative-norm waves and hence no particle creation [39] [60] • 

An interesting approach to this problem is to consider a discretized space whose 
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points are freely-falling, so that the spacetime is not strictly stationary and the lab 
frequency not necessarily conserved [ IU% , 109] . 

12.3 Final comments 

Despite some tantalizing experimental results - the realization of an analogue black 
hole in BEC [21], the detection of a signal from laser pulse filaments in nonlinear 
optical media [105] and the observation of the classical stimulated Hawking effect for 
surface water waves [22 02] - Hawking radiation remains stubbornly in the realms 
of theory. But the idea of its experimental realization is flourishing. Its concepts 
are constantly being applied to newer analogue systems. The goal is no longer a 
deeper understanding of gravity - and perhaps, given the surprising emergence of 
the “horizonless” regime, not so much about general black or white holes either. 
Instead, we are aiming to paint a picture of the quantum vacuum - of the content of 
nothing, of physics at its most fundamental [^J For example: the heuristic picture 
of particlc-antiparticle pairs arising from vacuum fluctuations and being torn apart 
by the event horizon might seem to give some insight into the gravitational black 
hole. But how is this picture to be reconciled with analogue white hole radiation, 
or the “horizonless” regime emerging from dispersion, where Hawking partners can 
be emitted in the same direction? Has our insight or our theory failed - and if so, 
what should replace it? For now, the vacuum and its contents might remain out of 
sight - but I hope I, and the many others in this held, can offer a tantalizing glimpse 
beyond the horizon! 


31 From this perspective, we can take encouragement from recent experiments on the related 
phenomenon of the dynamical Casimir effect (DCE), which predicts particle creation from the 
vacuum due to time-varying boundary conditions rather than a spatially-varying metric. Analogue 
systems in superconducting circuits of variable electrical length have been shown, both in theory 
|110l 1111] and experiment [11211113] . to produce analogue DCE radiation. 
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Figure 1: Radial WAVE TRAJECTORIES NEAR A BLACK HOLE HORIZON: Space- 
time diagram of light trajectories in the metric of Eq. ([5]), where the thick centre 
line represents the event horizon r = r$ and, in this near-horizon region, the flow 
velocity profile is approximately linear: V(r) ~ — c + a{r — r $). The dotted lines 
show trajectories of co-propagating waves (where the coordinate v, defined in the 
second of Eqs. (13), is constant), and these obey the second of Eqs. (|6]) ; these 
experience nothing unusual at the horizon, so that a wavepacket on the u-branch 
can pass through the horizon completely unhindered. In contrast, the solid lines 
show the trajectories of counter-propagating waves (where the coordinate u, from 
the first of Eqs. (13), is constant), and these obey the Erst of Eqs. © ; as these are 
propagated further back in time, they come to a standstill at the horizon. The u- 
branch is thus split into two pieces: Ur and ul (defined by Eq. (18)), its restrictions 
to the right- and left-hand regions, respectively. Both are right-moving with respect 
to the background fluid flow, but only the modes are able to propagate to the 
right in the lab frame, whereas the w^-modes are dragged to the left. Since the 
ul and Ur geodesics bunch together at the horizon, any wavepacket traced back 
in time will become arbitrarily thin, with arbitrarily short wavelength; this is the 
trans-Planckian problem. 
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Figure 2: BLACK HOLE HORIZON IN A RIVER: At the “Point of no return”, the flow 
speed of the water is exactly equal to the maximum speed of the fish. Any fish that 
travel further downstream than this point will be dragged inexorably towards the 
waterfall. [Illustration by Yan Nascimbene, from “Black Holes and the Information 
Paradox”, L. Susskind, Scientific American, April 1997] 
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Figure 3: WHITE HOLE HORIZON IN A RIVER: As in Figure [ 2 J the event horizon is 
the point at which the flow speed is exactly equal to the maximum speed of the fish. 
Here, however, the flow speed decreases in the direction of flow, so that the fish may 
swim arbitrarily close to the event horizon, but may never cross it. [Illustration by 
Peter Hoey, from “Test of Hawking’s Prediction on the Horizon with Mock ’White 
Hole’ ”, A. Cho, Science 319, 1321 (2008)] 
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Figure 4: DISPERSION PROFILE IN THE CO-MOVING FRAME: In the rest frame of 
the medium, the co-moving frequency obeys the simple relation uj 2 m = c 2 k 2 . This 
splits into four sectors: the re and rebranches are determined by the velocity in the 
co-moving frame, co cm /k = ±c (see Eq. (24)), and correspond to the diagonal lines 
ck and —ck, respectively; the positive- and negative-norm branches are determined 
by the sign of the co-moving frequency itself (see Eq. (29)), and therefore correspond 
to \ck\ (the thicker line on the upper-half of the plot) and — \ck\ (on the lower-half 
of the plot), respectively. The co-moving frequency is related to the conserved lab- 
frame frequency via the Doppler formula ca cm = u> — Vk, plotted here as dashed 
lines for the asymptotic velocities Vr and Vl and as dotted lines for the intervening 
velocities. The possible wavevector solutions occur at the intersections of the two 
curves, so that their evolution with changing velocity (i.e. with changing slope of 
the line oj — Vk) can be traced. Critical behaviour occurs as the velocity passes 
through —c, at which point the line u — Vk is parallel to the rebranch dispersion 
curve: it crosses the rebranch on the positive-norm side when V > —c, but on 
the negative-norm side when V < —c. As V —> —c, the wavevector on either the 
positive- or negative-norm side diverges - this is the trans-Planckian problem. 
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Figure 5: Dispersion profile IN the lab FRAME: Here is plotted the same 
information as in Fig. [4j but from the point of view of the lab frame. The conserved 
lab-frame frequency c o appears as a horizontal line. On the other hand, the dispersion 
profile of the medium in the lab frame is given by Vk ± \ck\, where the plus and 
minus signs refer to the positive- and negative-norm branches, respectively. The 
lab frame dispersion could also be written Vk ± ck , where this time the plus and 
minus signs refer to the u- and rebranches, respectively. As the flow velocity V is 
varied, the dispersion profile tilts, and the possible wavevector solutions - where 
the disperison profile is equal to w - vary accordingly. Again, we observe critical 
behaviour when V = —c, for at this point the u-branch of the dispersion lies exactly 
along the horizontal fc-axis, and is therefore parallel to the horizontal line at ui. For 

V > —c, it is the positive-norm w-branch which crosses oj, while for V < —c it is 
the negative-norm u-branch which does so. In contrast, the u-branch - which for 

V = —c is shown as a solid line within the continuum of V -dependent u-branches - 
experiences nothing of importance at V — —c. 
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Figure 6: SUBLUMINAL DISPERSION WITH SUPERSONIC FLOW: Here are plotted 
subluminal dispersion diagrams from the points of view of both the co-moving frame 
(upper panel) and the lab frame (lower panel). Positive- and negative-norm branches 
are shown as thick and thin curves, respectively, just as in Figs. [4] and |5} The velocity 
is assumed to be supersonic: V up < — Co, where Cq is the low-wavevector limit of 
the phase velocity, and therefore equal to the slope of the dispersion profile in the 
co-moving frame at k — 0. The line oj — V sup k, then, has a greater slope than the 
dispersion profile at k — 0, and since the slope of the dispersion profile decreases 
with increasing k, the two curves never cross for positive k. Viewed from the lab 
frame, the slope of the dispersion profile V sup k + c(k)k is negative at k = 0 because 
of the supersonic flow, and the gradual decrease of c(k) causes it to become more 
negative with increasing k, so that it can never equal c o for positive k. There are 
thus no solutions on the positive-norm u-branch, and only the negative-norm u -wave 
and positive-norm u-wave remain. 
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Figure 7: SUBLUMINAL DISPERSION WITH SUBSONIC FLOW: Here are plotted sub¬ 
luminal dispersion diagrams from the points of view of both the co-moving frame 
(upper panel) and the lab frame (lower panel). Positive- and negative-norm branches 
are shown as thick and thin curves, respectively, just as in Figs. ill and [6] The 
velocity is assumed to be subsonic: 14 u b > — Co, where Co is the low-wavevector limit 
of the phase velocity. As viewed from the co-moving frame, the slope of the line 
oo — KutA is therefore less than that of the dispersion profile at k = 0; but, since c(k) 
decreases with increasing k, the two curves become parallel at some point, and if u 
is equal to a; max the line u max — V su bk is actually tangent to the dispersion profile. 
As viewed from the lab frame, the dispersion profile Vk + c(k)k has positive slope at 
k — 0, but the slope decreases so that it eventually turns round; at this maximum 
point, it is equal to cn max . If uj < cu max , the two curves cross twice for positive k, 
giving rise to the ul- and w2-branches. The u-branch, then, does not extend to 
k — 0, but only to the point at which the (complex conjugate of the) u2-branch 
ends. On the other hand, if to > cu max , only the u- and u-waves exist. 
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Figure 8: DISPERSION WITH INHOMOGENEOUS V (SUPERSONIC TO SUPERSONIC): 
Here are plotted subluminal dispersion diagrams from the points of view of both 
the co-moving frame (upper panel) and the lab frame (lower panel). Positive- and 
negative-norm branches are shown as thick and thin curves, respectively, just as in 
Figs. [4j{7j The velocity is assumed to vary between two supersonic values, Ku P ,i and 
fsup, 2 - The entire spectrum is in the “high” frequency regime, where only the u- and 
v -waves - indicated by circles - exist. These vary smoothly with varying V, having 
well-defined values at all points in space; they do not experience a horizon. They do 
have opposite norm, however, so that they can still couple to form Hawking pairs. 
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Figure 9: DISPERSION WITH INHOMOGENEOUS V (SUBSONIC TO SUPERSONIC): 
Here are plotted subluminal dispersion diagrams from the points of view of both 
the co-moving frame (upper panel) and the lab frame (lower panel). Positive- and 
negative-norm branches are shown as thick and thin curves, respectively, just as in 
Figs. [4]{8j The velocity is assumed to vary between a subsonic value I4 u b > — Co and 
a supersonic value I4 up < — Cq. Therefore, of the critical frequencies, only cu maX! 2 is 
non-zero, and is labelled here simply as cj max . If lu > c u max , the situation is analogous 
to that described by Fig. [8j If co < cu max , we are in the mid-frequency regime 
where a group-velocity horizon exists. There, the ul- and w2-waves - indicated 
by inverted and upright triangles, respectively, in accordance with their lower and 
higher wavevectors - exist in the subsonic region, but they cannot exist in the 
supersonic region. Note their oppositely-directed group velocities - this is indicated 
by the slope of the dispersion profile relative to that of the line oj — Vk in the 
co-moving frame, or simply by the slope of the dispersion curve itself in the lab 
frame. If we track their evolution with varying V, we see that they vary towards 
each other, until, when V = V^ v h (indicated by a solid line among the continuum 
of curves), the straight line is tangent to the dispersion curve and the ul- and u2- 
wavevectors merge at the point indicated by a star. This point corresponds to the 
group-velocity horizon: as V varies further, the ul- and w2-waves cease to exist as 
real, propagating waves. Instead, they cross into each other, reversing their group 
velocity and moving back in the direction from which they originally came. They 
are thus not independent solutions, but degenerate into a single solution which 
“bounces” off the group-velocity horizon (see Fig. 12). This positive-norm solution, 
as well as the positive-norm r:-wave, can couple with the negative-norm w-wave to 
form Hawking pairs. 76 
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Figure 10: DISPERSION WITH INHOMOGENEOUS V (SUBSONIC TO SUBSONIC): Here 
are plotted subluminal dispersion diagrams from the points of view of both the 
co-moving frame (upper panel) and the lab frame (lower panel). Positive- and 
negative-norm branches are shown as thick and thin curves, respectively, just as in 
Figs. [4]j9j The velocity is assumed to vary between two subsonic values, I4ub,i and 
Vjjubg. The critical frequencies o w-r i and are therefore both non-zero, so that 

all three velocity regimes described in §8.2| exist. It is easiest to pick them out in 
the lab frame diagram: the critical frequencies are there simply the maxima of the 
dispersion curves corresponding to the asymptotic velocities. If cu > co max , 2 , then 
only the u- and u-waves exist, so that the situation is analogous to that of Fig. [8j If 
w m ax,i < (jJ < cu maXj 2 , then the ul- and w2-waves only exist in one of the asymptotic 
regions, and the situation is similar to that described in Fig. 0 « lu < u; maXj i, 
then the ul- and w2-waves exist in both asymptotic regions, varying smoothly as 
V is varied. They do not degenerate into a single solution, so that there are four 
independent solutions in total. The ul-, u2- and u-waves all have positive norm, 
and can couple with the negative-norm u -wave to form Hawking pairs. 
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Figure 11: wl-OUT MODE IN BLACK HOLE CONFIGURATION: Here is illustrated the 
wl-out mode in the presence of a black hole horizon with subluminal dispersion. The 
left panel shows a space-time diagram of the trajectories of the various wavepackets, 
which are labelled with the same symbols used in the dispersion diagram of Fig. [9j 
on which the evolution can be traced. Included as an inset is a close-up of the group- 
velocity horizon region, where the incoming w2-wave crosses onto the wl-branch at 
the point marked by a star, “bouncing” off the group-velocity horizon. Note that 
the simple geometric picture provided by the dispersion diagram cannot describe 
coupling into the u- and w-branches. (The trajectories of the u- and r:-waves are 
almost equal near the group-velocity horizon, making them difficult to distinguish). 
In the right panel are shown the results of numerical wavepacket propagation at three 
representative times: the asymptotic past, the interaction with the near-horizon 
region, and the asymptotic future. This was calculated by specifying the outgoing 
ul-wave and propagating it backwards in time. (For these calculations was used 
the dimensionless velocity and dispersion profiles of §|ll.l C 2 (K ) = 1 — K 2 and 


U(X ) = {Ur + Ul)/2 + (Ur — Ul)/2 tanh(aX), with Ur = —0.5, Ul = —1.5 and 
a = 1.) While the outgoing ul-wave has positive norm, the u-wave is the only 
ingoing wave with negative-norm; so, according to 19.2 the norm of the u -wave 


relative to the wl-wave gives the creation rate of the outgoing wl-wave. 
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Figure 12: w-IN MODE IN HORIZONLESS CONFIGURATION: Here is plotted the u -in 
mode in the absence of a group-velocity horizon with subluminal dispersion. The 
flow is entirely subsonic, and the frequency is in the low-horizon regime c o < cu maX) i, 
so that the ul- and w2-waves exist everywhere. The situation is thus similar to that 
described in Fig. [lOl although the simple geometric picture provided by wavevec- 


tor evolution on the dispersion diagram cannot account for coupling between the 
various waves. The left panel shows a space-time diagram of the trajectories of 
the various wavepackets, which are labelled with the same symbols used for the 
various wavevector solutions in Fig. 10 The right panel shows the results of nu¬ 


merical wavepacket propagation at three representative times: the asymptotic past, 
the interaction with the inhomogeneity in V, and the future. Note that the u2- 
and w-waves have almost the same group velocity, so that they lie on top of each 
other in the future and are difficult to distinguish. (For these calculations was used 
the dimensionless velocity and dispersion profiles of §11.1 C 2 (K ) = 1 — K 2 and 


U(X) = (U R + U L )/2 + (U R - U l )/2 tanh(aX), with U R = -0.4, U L = -0.8 and 
a — 1.) The u-wave is the only one with negative norm; therefore, according to 
the norms of the outgoing ul-, u2- and r^-waves relative to that of the ingoing 


).2 


u -wave give their respective creation rates. 
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Figure 13: INTEGRATION OF STATIONARY SOLUTION: In one asymptotic region, 
the solution is chosen to be a particular plane wave. Eq. (72) is then integrated 
through to the other asymptotic region, where the solution is a linear combination 
of plane waves (73). This is done for all possible plane wave solutions in the initial 
region. Thus we End the transfer matrix T appearing in Eqs. (74) and (75). 
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Figure 14: PHASE INTEGRAL IN MOMENTUM SPACE: Using the dispersion and 
velocity profiles of 1 11.1| (plotted in Fig. 15), here the colour shading shows the 
imaginary part of the position y(fc) - the darkest shade represents —7r/2a, the 
lightest +7r/2a - where the wavenumber k is taken as a complex variable. The 
dashed blue curves show where \{k) is real, and split into two disconnected regions: 
the cross-bow shape on the right contains (on the real axis) the positive-norm ul- 
and w2-branches, as well as the complex conjugate wavevectors they connect to 
beyond the group-velocity horizon; the line on the left contains the negative-norm 
w-wavevectors. The extremities of these curves correspond to the solutions in the 
asymptotic constant-velocity regions; those on the left are shown in black, those 
on the right in red, and the arrows show the directions of steepest descent through 
these saddle points. When solving the wave equation, the boundary condition is that, 
beyond the horizon, only the exponentially decreasing solution is allowed; therefore, 
for large negative x, the integration contour must be able to pass through the saddle 
point near k^ c in its direction of steepest descent - this contour is shown in black, 
and does not contain any contribution from the other left-hand (colour black) saddle 
points. As x varies continuously into the right-hand asymptotic constant-velocity 
region, this contour is continuously deformable to the red one, which is able to pass 
through all three of the right-hand (colour red) saddle points in their directions 
of steepest descent. The relative phases between the various plane waves is given 
by the integral of y(/c) between their corresponding wavevectors, along the path 
selected by the integration contour. When comparing the positive- and negative- 
norm components, this integral traverses a branch cut of y(/c) and hence contains 
an imaginary part, which can be isolated from the real part by taking the closed 
contour integral around the entire branc^qcut. This gives their relative amplitudes 
- and the Hawking temperature - as described by Eqs. (95) and (96). 



















Figure 15: SPECIFYING DISPERSION AND VELOCITY PROFILES: 

(a) The solid line plots c(k)k , with c(k) as in Eq. (98). The dotted line plots the 
corresponding dispersionless curve, cok. 

( b ) The velocity profile, given by Eq. (99), is monotonic and approaches asymp¬ 
totically constant values. The parameter a is inversely related to the length of the 
transition region. 
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Figure 16: TEMPERATURE - VARYING U L (LOW -a REGIME): The (normalized) 
Hawking temperature, for the dispersion profile (102) and the velocity profile (103), 
is shown as a function of frequency for various values of Ul, while Ur is fixed at 
—0.5 and a is fixed at 0.1. The solid curves correspond to the numerically calcu¬ 
lated spectra of the wl-wave, while the thin dashed curves show the spectra for 
the u2-wave, with a maximum frequency (indicated by vertical dotted lines) that 
vanishes when \U > 1. The thick dashed curves show the analytic prediction of 
Eq. (105), derived from the phase-integral method described in f 10.3 this is applied 
only when a group-velocity horizon is present, and so only for the wl-wave spectra. 
The agreement is seen to be very good for a = 0.1 and for Ul < —0.7. (It is not 
so good for Ul — —0.51, where the group-velocity horizon exists only over a very 
narrow frequency range.) 
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Figure 17: TEMPERATURE - VARYING U L (DISCONTINUOUS LIMIT): The (normal¬ 
ized) Hawking temperature, for the dispersion ( |102 ), in a step-discontinuous flow is 
shown as a function of frequency for various values of Ul, with Ur fixed at —0.5. 
The solid curves correspond to the wl-wave, the dashed curves to the w2-wave (cou¬ 
pling into which only occurs in the purely subsonic regime where \Ul\ < 1), and 
we have also included dotted lines showing the temperature of the u-wave, which is 
easily calculated in the step-discontinuous limit. (The highest u-curve corresponds 
to Ul = —1.5 and the lowest to Ul = —0.51, with monotonic variation in between.) 


The most significant difference with respect to the low-steepness regime of Fig. 16 


is in the temperature values, which are here (for ul and u2) an order of magnitude 
greater. The shapes of the spectra are essentially the same as in Fig. 16, except 
that the fall-off on the approach to fl maXi 2 begins to take hold at lower frequencies, 
and the connection between the regimes in which a group-velocity horizon does and 
does not exist is a more clearly visible feature in the spectra of the wl-wave. The 
r?-wave spectra do not go to zero at f2 max because the u-v coupling (between positive 
and negative norm) exists even above this value; this can be seen in Figs. [6] and [9j 


84 



























Figure 18: TEMPERATURE - VARYING Ur (low -a regime): The (normalized) 
Hawking temperature, for the dispersion profile ( 102[ ) and velocity profile (103), 
is shown for various values of Ur while Ur is fixed at —1.5 and a is fixed at 0.1. 
Since Ur < —1, there is no “horizonless” regime and hence no emission of u2- 
waves. The solid curves show the numerically calculated wl-wave spectra, while the 
dashed curves correspond to the analytical prediction of Eq. (105). The maximum 
frequency (shown as dotted vertical lines) decreases with Ur, eventually vanishing 
in the limit Ur —> — 1. 
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Figure 19: TEMPERATURE - VARYING Ur (DISCONTINUOUS LIMIT): The (normal¬ 
ized) Hawking temperature, for the dispersion profile (102) and a step-discontinuous 
flow, is shown as a function of frequency for various values of Ur while Ur is fixed 
at —1.5. Solid lines show spectra for the wl-wave, and dotted lines for the u-wave; 
since Ur < —1, there is no “horizonless” regime and no emission of u2- waves. Ver¬ 
tical dotted lines show the maximum frequencies at which the ul-u coupling ceases 
and the spectrum vanishes; since the u-v coupling can occur at all frequencies, the 
r?-wave spectra do not experience such a cut-off. As in the case of varying Ur (see 
Figs. 16 and m. we find in comparison with the low-steepness regime of Fig. [18 


that the temperatures are greater by an order of magnitude and the fall-off is more 
noticeable at lower frequencies. 
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Figure 20: LOW-FREQUENCY TEMPERATURE - VARYING a: Here is shown the 
(normalized) low-frequency temperature with increasing steepness a, for various 
values of Ul while Ur is fixed at —0.5. The parameter a is normalized with respect to 
its value at which the linear Hawking prediction is exactly equal to the temperature 
in a step-discontinuous flow; they are equal at the dotted line, and the value of a 
in relation to this gives a good indication of the point of transition between the two 
regimes. Dashed lines correspond to the linear Hawking prediction for low a and 
the limiting step-function temperature for high a. 
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